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ANALYSIS  AND  MINIMUM  WEIGHT  DESIGN 
IN  SANDWICH  BEAMS  WITH  RIGID  FOAM  CORE 


I.  Introduction 


A  sandwich  beam  is  shown  in  Figure  1.  It  has  a  span  i  and  a 

width  b.  The  face  and  core  thicknesses  are  t  and  c  respectively 

and  the  core  dens i tv  is  p  .  Tne  face  material  has  a  density  p  .  a 

c  ^  i 

Young's  modulus  E^  and  a  yield  strength  o  The  cote  has  a 

Young's  modulus  E  ,  a  shear  modulus  G  ,  a  uniaxial  yield  strength 

*  • 

a  and  a  shear  strength  r  .  The  solid  from  which  the  core  is 

C  C 

foamed  has  a  density  p  ,  a  Young's  modulus  E  and  a  yield  strength 
c 

yi 

In  order  to  simplify  the  analysis  of  sandwich  beams,  the 
following  assumptions  are  made  in  this  program  : 

(1)  Both  the  face  and  core  materials  are  isotropic. 

(2)  The  sandwich  beam  has  an  antiplane  core  and  thin  faces. 

(3)  The  beam  is  narrow,  so  that  normal  stresses  in  the  y 
direction  can  be  neglected. 

(4)  Ordinary  beam  theory  is  valid  for  the  sandwich  beam. 

The  foam  core  property-density  relationships  are 


E  -  C  (  p  /  p  )"  E 

c  3  e  •  * 


r  -  C  (  p  /  p  )  o 

e  *  e  »  y» 


G  -  C  {  p  /  p  )  E 

c  g  c  s  • 


where  C  ,  C  ,  C  ,  A,  B  and  G  are  material  properties, 

3  *  t 


1 


For  a  sandwich  bean,  the  maximun  moment  is  given  by 
Pi 

M - 

C 

1 

The  maximum  shear  force  in  the  bean  is 
P 

C 

2 

The  maximum  deflection  due  to  bending  and  shearing  stresses 
can  be  expressed  as 

t  -  L  +  L 

b  t 

pi^  Pi 


CD  C  A  G 

i  6  •  c 

where  C  ,  C  ,  C  and  C  are  constants  related  to  the  loadinr 
i  2  i  6  ^ 

geometry. 

bc^  btd^  bc^ 

D  -  E  -  +  E  -  +  E  - 

'  6  '2  '12 


bd 


2 


A 


c 


d  -  c  +  c 

Figure  2  shows  the  load  constants  for  six  different  loading 
configurations  which  are  available  in  this  program. 

The  program  has  four  parts  : 

(1)  Analysis  of  sandwich  beam  deflections  and  stresses. 

(2)  Minimum  weight  design  of  a  sandwich  beam  for  stiffness. 

(3)  Minimum  weight  design  of  a  sandwich  beam  for  strength. 

(4)  Minimum  weight  design  of  a  sandwich  beam  for  both 
stiffness  and  strength. 
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II,  Analysis  of  Sandwich  Bean 


The  critical  normal  and  shear  stresses  in  the  face  material 


a - E 

*  D  ' 


(  c/2^z^/z  :  -h/2<2<-C/i  ) 


- - E  ( 

f  t  ‘ 

D 


((c/2)+C)*  -  Z* 


The  critical  normal  and  shear  stresses  in  the  core  are 


a  -  -  E 

'  D  ' 


(  -C/2<2<C/2  ) 


V  td 


E  2 

c  C 

-  (  -  •  2^  )  ] 

2  4 


The  maximum  deflection  in  the  beam  is 


CD  C  A  G 

i  6  •  c 


Note  that  the  location  of  applied  force  can  be  arbitrary’  in 
the  analysis  of  sandwich  beam. 

The  failure  equations  for  the  sandwich  beam  are  : 

Face  yielding 


P  -Co  bct/l 
tr  I  yi 


Face  wrinkling  : 


2/3  1  ^  '3  2A/3 

-  0.57C  C,  V  -E  “’(p  /p  )^'hct/l 

13Cs  ct 


3 


Core  shearing  : 


P 

ca 


C  {p  /p  )“ff  he 

4  e  •  7  ■ 


CAfi  /P  )*E 
3  6a 


a 


{  [ - 

2C^(t/i)E^ 


1  /  2 

) 


This  program  vill  make  a  comparison  of  P  .  P  ,  P  and 

^  fy  £w  e» 

pick  one  critical  loading,  then  check  whether  the  stress  state  is 
safe  or  not. 


Ill,  Minimum  Weight  Design  for  Stiffness  in  Sandwich  Beams 


The  maximum  deflection  of  the  sandwich  beam  is 

Pi’  Pi 

^  .  -  + - 


CD  C  A  G 

i  6  •  c 


Solving  for  the  density  of  the  core,  we  obtain 


ICE 
^  j _ s  _ f 

C 


Pitc 


l/G 


C  C  E  (  AC  bze^E  -  2Pi’  ) 

»  6  •  5  t 


The  weight  of  the  sandwich  beam  is 


V  -  2p  biz  +  p  blc 

i  e 


ICE  1  /  G 

-  2p  bci  +  ( - ’  , 

‘  C  C  E  ‘ 

t  6  s 


[  AC  bcc^E  -  2Pi’  ] 

^  a  »  * 


-l/G 


Letting  8'^/8t  -  0  and  SV/Sc  -  0  and  solving,  we  find  that 


G  4  :/G  2+2G  i-:/G 

c  .  -  (  4[ -  ( - )  ]( - ) 

G-1  G-1  G-1 


iC  PE  , 

5  f  -1/G  1  1-1/C 

( - )  T£\ - ) 


C  C  Z  AC  bE 

(  6  •  5  i 


C/(3G-1) 
)  ■ 


2(l+G)Pi' 


(G-l)AC  E  be  ^ 

i  f  opt 


iC  PE  t  c 

5  f  opt  0  p 1 


l/C 


(P  )  -  P  [- 


‘  ‘  C  C  E  (AC  E  be  c  ‘  -2Pi’) 

ge«  optopt 


Also,  we  can  calculate  the  r.cndimensional  parameter  which 
gives  a  measure  of  the  relative  face  and  core  stiffnesses. 


C  (p  /p  )'•  E  c  ^ 

t  c  (opts  opt 


opt  1/2 


B  - 


(1+3-——)  ) 


opt 


2E 


opt 


opt 


Allen  suggests  that  B  should  be  greater  than  20  tr  ensure 
that  shear  lag  does  not  occur. 


IV.  Minimum  Weight  Design  for  Strength  in  Sandwich  Beams 


There  are  three  possible  failure  modes  considered  in  this 
program.  For  each  case,  the  critical  loading  can  be  calculated 
Face  yielding  : 


5 


P  -  C  cr  bct/i 
ty  1  yf  ' 


Face  wrinkling  ; 


2A/3 


P  -  0.57C  C  '  E  E  (  p  /P  )  bct/i 

Iw  1  3  £  «  e  • 


Core  shearing  ; 

Assuming  a  «  r 


P  “  C  C  <7  (  P  /p  }  be 

et  2  *  7>  c  a 


The  transition  equation  between  two  failure  modes  can  be 
found  by  equating  the  corresponding  applied  forces. 

Face  yielding-Face  wrinkling  : 

P  -  P 

£y  £« 


y£ 


3/2A 


•>  P  /P  -  ( 


'  •  0.S7C  ‘"e  ‘'’e  = 

3  £  t 


Face  yielding-Core  shearing 


P  -  P 

£y  ca 


c  c 

2  * 


P  c 

B  y  a 


C/i  - 


(  -  ) 


1  y£ 

Face  wrinkling -Core  shearing  : 


P  -  P 

£w  ca 


c  c 

2  « 


— >  C/i  - 


0.57C  c/'’  P 

13  a 


B-2A/3 

(  -  ) 


y  a 


E 


A  typical  failure  mode  map  is  shown  in  Figure  3 . 
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The  optiiEuni  design  lies  on  the  transition  line  betveen  two 
failure  modes.  There  are  four  possibilities . 

(a)  Face  yielding-Face  wrinkling  fail  simultaneously 


P  -Co  bct/i 
iy  lyt 


Pi 


— >  c  - 


Co  be 
1  yt 


P  -  0.57C  C  /p  )^"bct/2 

Jw  1  3  f  »  c  f 


2A/3. 


y  i 


3/2A 


- - ( 

P 


0.57C  ^^'E  ^'“E 


3  !  i 


Therefore,  the  weight  function 

Pi 


y f  3/2A 


W  -  2^  bi( 


Co  be 

1  yt 


)  +  ( 


)  P  bei 

f 


,  „  c-»/-  l/3_  2/3 

where  Q-0.57C  E  E 

3  r  » 

Letting  dV/dc  -  0 ,  we  find 


2p  Fia 


3  /  2A 


1/2 


C  -  [ 

opt 


-  ,  .(1+3/2A), 

Co  (o  )  b 

It  y  f 


Pi 


Opt 


Co  be 

1  yt  opt 


t  opt  t  Q  2/3g  l/3g  2/3 

3  {  1 


3/2A 

)  -  constant 


From  the  failure  mode  map,  it  is  realized  that  there  is  a 


/ 


restriction  in  c/i  for  the  co-occurrence  of  face  yielding  and  face 
wrinkling  failures. 


opt 


C  C 

t  2  * 


a 

B  r  I 


C  a  be  1  C 

1  y  r  opt  1 


{p  IP  )  , 

c  •  opt 


yf 


-  C  *  2p 

P  2  *  t 

->  -  <  -  a  ^  -  ( 

bi  C  '''  p  a  _  a 


a  (6B-3)/2A 

yf 


1 


•  yf 


(b)  Face  wrinkling-core  shearing  fail  simultaneously 


P  -  0.57C  C  /p  )^'^bct/l 

fw  I  3  f  (  c  « 


Pi 


— >  c  - 


0.57C  C  *'’e  ^'’e  ‘'*(p  /p  )^*'''bc 

13  1  »  e'  • 


P  -CC(p/p)a  be 

c*  2  *  e  •  y« 


— >  p  ^  p  { 

e  f 


P  l/B 

-  ) 


C  C  <7  be 

2  *  y  * 


Therefore,  the  weight  function 

.2 


Pi‘ 


W  -  2p. 


C  a(p  /p  )*”'  c 

1  e  • 


-  +  p  i(- 

2A/3 


P  i/B  1-1/B 

- )  (be) 


C  C  e 

2  ♦  y» 


Letting  aw/3c  -  0.  It  is  found 


p  (3-3B)C  ob 

1  P  (2A+3)/3B  3B/(2A-6B+3) 


C  - 

opt 


2p^(2A-3B)Pi  C^C^o^b 


■)  ] 


8 


(p  )  ~  p  ( 

e  opt  1 


C  C  a  be 

2  *  yt  opt 


C  a(p  /p 

1  c  •  opt  opt 


Similarly,  from  the  failure  mode  map,  there  is  a  restriction 
in  p  /p  for  face  wrinkling-core  shearing  failure  mode. 


{p  /p  )  p  /p  -  ( 

c  t  opt  e  ( 


yt  3/2A 


Therefore , 


2p  (2A-3B)C  ^ 

i  2  *  y» 


/j^(3-3B)C^o 


a  (6B-2A-3)/2A 

yi 


(c)  Face  yielding-core  shearing  fail  simultaneously 


P  •  C  e  bcc/l 
ty  1  yt 


— >  c  - 


C  O'  be 
1  yf 


P  -  C  C  (p  /p  )*a  be 

ei  2  4  c  •  ys 


•>  p  ^  p  (, 

c  ■ 


C  C  o  be 

2  «  7* 


Thus,  the  weight  function  is  : 


9 


V  - 


2pYi‘ 


C  a  c 
1  yf 


C  C  a 

2  *  y  ■ 


1/B  1-1/B 

)  (be)  p  i 


Letting  cV/Sc  -  0.  L'e  can  find 


2p^BPi 


C  C  ff  b 

2  *  yi  1/B  B/(2E-1) 


c  -[ 
opt 


p  (B-l)bC  a 
«  1  y  f 


)  ] 


Pi 


t  - 

opt 


Co  be 

1  yi  opt 


1/t 


(P  )  -  P  ( 

c  opt  s 


C  C  a  be 

2  *  y»  opt 


In  order  for  face  yielding  and  core  shearing  to  occur  at  the 

same  time  there  is  a  limitation  in  p  /p  according  to  the  failure 

c  • 

mode  map.  That  is 


yi  3/2A 

(p  /P  )  ^  P  /P  -  (  -  ) 

c  t  opt  c  t 


Therefore , 


„  2p  BC  ‘a  ^ 
P  I  2  *  yi 


yf  (6S-3)/2A 
(  -  ) 


bi  p  (B-l)C  a  ^ 
•  1  yf 


(d)  Face  yielding- face  wrinkling  -core  shearing  fail 


simultaneously 


P  -Co  bec/i 
fy  1  yf  ' 

P^  -  C  Q(p  /p  )“^^bcc/i 

f«  1  C  1 
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P  -  C  C  a  {p  /P  ybc 

es  2  4  yt  e  • 


The  three  equations  in  three  unknowns  give  : 


P 


c 

opt 


C  C  <7  b 

2  4  jr  t 


a 

yt  -3B/2A 

-  )  -  constant 

a 


c  ' 

opt 


Pi 


C  a  be 

1  yt  opt 


constant 


a 

y  i  3/2A 

(p  )  “  p  (  -  )  "  constant 


In  running  this  program,  the  P/bi  requirement  is  checked 


first.  The  optimal  design  is  evaluated  for  each  case  for  which  the 


condition  is  met.  Among  the  above  possible  failure  modes,  the  one 


with  the  minimum  weight  is  selected  as  the  overall  optimum  design. 


V.  Minimum  Weight  Design  for  Both  Stiffness  and  Strength 


In  practice,  the  minimum  weight  design  should  take  both 
stiffness  and  strength  requirements  into  account.  In  this 
program,  there  are  three  kinds  of  failure  mode  considered  for 
sandwich  beams. 

(a)  Stiffness  and  face  yielding  failure 


P  -  C  bec/i 
ty  I  yt 
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Pi 

— >  c  -  - 

C  a  be 
1  yt 


Pi 


C  A  G 

6  •  c 


C  E  Pitc 

5  t 


— >  P  -  [ 


C  C  E  (C  Abec  E  -2Pi  ) 

t  6  •  5  f 


1/G 
1  P 


Therefore,  the  weiphe  function  becomes 


W  -  2p  bti  + 


AC  C  E  C 

(  6  ■  1 


(G-1 )  (G-  1  ) 

O 

yi 


6  ,2S 
i 


^(G-i) 


2C  C  E  C 


(6(1 


yf 


C  E  o^«2G-n 

it  * 


-1/6 

1 


Since  the  weight  of  the  sandwich  beam  is  a  real  number, 
the  value  in  the  bracket  should  be  larger  than  zero.  From  this 
observation  it  is  concluded  that 

AC  E  P 

5  f 

0  <  c  <  - 

2C  *bf 
1  yt 

There  is  a  constraint  for  face  yielding  failure. 


<P  /P  )  , 

c  s  opt 


p  /p 

c  t 


( 


O 

yi 


a 


3/2A 

) 


Within  the  range,  a  minimum  weight  can  be  found  by 
incrementally  increasing  t.  The  corresponding  face  and  core 
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thicknesses  and  foam  density  are  the  optimal  design  values. 


(b)  Stiffness  and  face  wrinkling  failure 


p  - 

fw 

c  q(p  /p 

I  e 

)^^hct/£ 

1 

— > 

pf 

.  .  .-2A/3 

(  P  /P  ) 

c  s 

C  abc 

1 

Pi^ 

Pi 

A  “ 

C  D 

j 

C  A  G 

6  •  c 

2i‘c  a 

A 

Pi 

,  .2A/3  . 

■  (P  /P  )  + 

‘  ‘  C  C  E  bA 

c  • 

C  AE 

i  K  6  s 


Tnus ,  the  weight  function  can  be  expressed  as 


2fi  bPi^  C  C  E  C  AE  (p  /p 

I  gCtSfct 

C  a  ^  2i‘C  C  E  oC  b*C  E  Pi(p  /p  )  ‘ 

1  (StlSf  cs 


2i^C  C  E  C  Qb  +  C  E  Pi(p  /p  )‘'2A'3C)/3 

(6il  if  c( 

C  C  E  C  AE  (p  /p 

(  6  s  i  f  C  s 


From  the  failure  mode  map,  it  is  known  that  there  is  a 
limitation  in  p  /p  for  occurrence  of  face  wrinkling. 

C  S 

a 

0  <  P  /p  <  (  — ^ 

c  « 

a 

There  is  a  constraint  for  face  wrinkling  failure. 


t  E-(2A/3) 

-  <  -  (p  /p  ) 

.  -  e  • 
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The  value  which  corresponds  to  the  minimuni  weight  can  be 
found  using  the  sane  argument  as  that  for  the  constraints  of 
stiffness  and  the  face  yielding  failure  mode. 

(c)  Stiffness,  face  yielding  and  face  wrinkling 


P  -Co  bct/i 
fy  i  yt  ' 


P,  -  C  Q(p  /p 

1«  1  e  ( 


Pi 

t  -  -  +  - 

CD  C  A  G 

i  e  •  c 

The  three  constraints  in  three  variables  can  be  solved  for 


yl  *3G/2A 
—  ) 


Pi 

Co  be 

1  yi  opt 

In  running  this  program,  the  three  failure  modes  considered 
above  will  be  compared  with  each  other  to  determine  the  overall 
minimum  weight  design. 

It  is  noted  that  the  stiffness -core  shearing  failure  mede  is 
not  included  in  this  program  because  it  ends  up  with  an 
impractically  large  core  thickness. 


yi 


3/2A 


iP  /P  )  -  ( 

e  t  opt 


c 

opt 


2i  C  c 
1  yi 


C  AE 

5  r 


Pi 


C  C  E  bb 

S  6  ■ 
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Figure  2.  The  load  constants  for  different  geometry. 


CORE  RELATIVE  DENSITY. /7/p^  (-) 


I0‘'‘  lO"’  lo'^  10“' 


FACE  THICKNESS /SPAN  LENGTH,  t/l  (-) 


Figure  3.  A  failure  node  nap  for  a  sandwich  bean  loaded 


in  three-point  bending  with  strength  contours 
superinposed  on  the  nap  (  b-c-25  mn  ) . 
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ANALYSIS  AND  MINIMUM  WEIGHT  DESIGN 


IN  SANDWICH  PLATES  WITH  RIGID  FOAM  CORE 


I.  Introduction 


There  are  three  kinds  of  sandwich  plate  considered  in  this 
program  : 

1.  A  simply  supported  circular  sandwich  plate  under  a 
distributed  load  over  a  central  circular  portion  of  its 
area . 

2.  A  clamped  circular  sandwich  plate  under  a  distributed 
load  over  a  central  circular  portion  of  its  area. 

3.  A  simply  supported  rectangular  sandwich  plate  under  a 
distributed  load  over  its  entire  area. 

Figure  4  shows  a  circular  sandwich  plate  which  spans  a  radius 
r  and  carries  a  distributed  load  q  over  a  central  circular  portion 
of  its  area  of  radius  a.  The  thickness  of  each  face  is  c,  while 
that  of  the  foam  core  is  c.  The  faces  have  a  density  f>^,  a  Young's 
modulus  E  ,  Poisson's  ratio  w  and  a  yield  strength  a  .  Tne  solid 

polymer  of  which  the  core  is  made  has  a  density  p  .  a  Young's 

modulus  E  and  a  yield  strength  a  .  The  density  of  the  foamed 

core  is  p  ,  its  Young's  modulus  is  E  .  its  shear  modulus  is  G  and 

e  e  e 

its  shear  strength  is  . 

A  rectangular  sandwich  plate  which  has  a  length  b  and  width 
a  (  note  that  b  "i.  a  ) ,  is  shown  in  Figure  5. 

There  are  three  assumptions  made  in  this  program  : 

(1)  Both  the  face  and  core  materials  are  isotropic. 
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(2)  The  sandwich  plate  has  an  antiplane  core  and  thin  faces. 

(3)  Stresses  in  the  faces  and  core  in  the  z-direction  are  of 
no  importance  and  neglected. 

The  foam  core  property-density  relationships  are 


E 


C 

3 


{  P  /  P  )*  E 

C  S  ft 


r  -  C  (  p  Ip  )  a 

c  *  c  s  y> 


G  -  C  (  p  /  ^  )“  E 

C  g  C  ft  ft 


where  C  ,  C  ,  C  ,  A,  B  and  G  are  material  properties. 

For  a  sandwich  plate,  the  maximum  deflection  can  be  expressed 


as 


qa. 


qa 


t  - 


16D 


g  + 


AS 


The  maximum  stresses  in  the  face  material  are 


qa 

<  or  ) - (  ) 

cc 


qa 

a  (  or  c  ) - (  g^  +  ) 

cc 


qa 

r  (  or  r  ) - (  1  *  *'-  )  g. 

^  rt  cc  '  ^ 


The  maximum  stresses  in  the  core  are 


f 

XZ 


r 

yt 


(  or  r 

rt 


(  or  r 

tx 


) 

) 
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are 


constants 


where  g^.  g^.  g^.  g^.  g^,  g^,  g^,  g^  ^6  g^ 

related  to  the  loading  geometry  of  the  sandwich  plate. 

E  c 

f 

D - 

2  (  1  -  ) 

S  -  c  G 

c 

The  loading  geometry  constants  for  three  different 
configurations  which  are  available  in  this  program  is  given 
in  Appendix. 

The  program  has  four  parts  : 

(1)  Analysis  of  sandwich  plate  deflections  and  stresses. 

(2)  Minimum  weight  design  of  a  sandwich  plate  for  stiffness. 

(3)  Minimum  weight  design  of  a  sandwich  plate  for  strength. 

(4)  Minimum  weight  design  of  a  sandwich  plate  for  both 
stiffness  and  strength. 


II.  Analysis  of  Sandwich  Plate 


The  maximum  stresses  in  the  face  material  are 

2 

qa 

(  or  C7^  )  -  -  (  g^  +  ) 

*  *  cc 

2 

<7« 

<7^  (  or  )  -  -  (  g^  +  ) 

2 

qa 

> - (  1  -  )  s, 

CC 

The  maximum  stresses  in  the  core  are 
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r  (  or  r  ) - g 

X2  r*  8 

c 


T  (  or  r  ) - g 

yi  ti  °i 

C 


The  maximum  deflection  in  the  plate  is 


-  g,  +  - 

16D  4S 


The  failure  equations  for  the  sandwich  plate  are 


Face  yielding 


c  c  t 
yt 

c  .  - 

'ty  2  ,  . 

^  (  F3  +  ) 


Face  wrinkling  : 

Plantema  pointed  out  that  the  superposition  of  the 
perpendicular  normal  stresses  in  the  plane  of  the  face  does  not 
affect  the  wrinkling  stress.  Therefore, 


a  c  t  {p  /p 

c  ■ 

C  -  - 

a  <  S3  +  ) 


,  -  2/3  1/3  2/3 

where  o  -  0.57C  E  E 

3  f  » 

Core  shearing  : 


Assuming  a  «  r  and  the  core  material  is 


isotropic.  Thus, 


(P/P.) 


This  program  will  make  a  comparison  of  q^^  and  pick 
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one  critical  loading,  then  check  whether  the  stress  state 
or  not. 


III.  Minimum  Weight  Design  for  Stiffness  in  Sandwich  Plates 


The  maximum  deflection  of  the  sandwich  plate  is 

2 


ca 


qa 


A  - 


16D 


4S 


Since 


G  -  C  (fi  /fi  r  E 

C  (  c  *  s 


Solving  for  the  density  of  the  core,  it  is  found 


2  .  * /I  2, 

a  A  a  (l-i/^  )g^ 


P  “  P  { 

e  t 


(  -  - 


-1  :/c 

]  ) 


AC  E  c  o  BE  tc 


The  weight  of  the  sandwich  place  is 

V  -  2p  abc  +  p  abc  (  or  W  -  2p  nr^c  +  p  rr^c 

t  c  {  c 

Letting  dV/8c  -  0  and  dV/Bc  -  0  and  solving,  we  find 


2  ■  f  qa  2  X/(3G'l) 


(p  )  -  p  I 

c  opt  • 


(  —  )  ) 


256C  G(1-k  )gpZ  A 

»  t  I  « 


<-5  ,  vC*  1  ,,  2.C  c  G_  1 

a  °  ^  ^  qa  g-i  — 

c  -  -{A®*^C  -  (— ) 

‘*^*'2  ‘  ,  ,.2G  G_  G  A 

(G-1)  Ej 


safe 


) 
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c  -  4a 

opt 


,_2  , 2.G-1  C- J  2_  :-G 

(G  -1)  g^ 


2l3G*lc  Z^2G^  2 


qs  1/(36'1) 

(p  /p  )  (  -  ) 


IV.  Miniaun  Weight  Design  for  Strength  in  Sandwich  Plates 


Similar  to  those  of  the  sandwich  beams,  there  are  three 

possible  failure  inodes  considered  in  this  program. 

Face  yielding  : 

o  c  t 
ri 


^iy  “  2  . 

^  ) 

Face  wrinkling  : 

a  c  t  {p  /p  ) 

e  • 


^fw  "  2 


Core  shearing  : 


Co  c  {p  /p  ) 

♦  y  1  c  » 


- 


«  s. 


The  transition  equation  between  two  failure  modes  can  be 
obtained  by  equating  the  corresponding  applied  distributed  forces. 
Face  yielding- Face  wrinkling  : 
o  -  q 


— >  P  /p 

c  t 


O 

y  1 


Q 


3/2A 

) 
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Face  yielding-Core  shearing  ; 


£y  e» 


->  C/a  - 


o  g 


Face  wrinkling-Core  shearing  ; 


<7,  -  <7 

!<•  e» 


B  -  2A  /  3 


The  failure  mode  map  is  the  same  as  that  of  sandwich  beams 
except  that  the  coordinate  t/i  replaced  by  c/a. 

The  optimum  design  lies  on  the  transition  line  between  two 
failure  modes.  There  are  four  possibilites . 

(a)  Face  yielding-Face  wrinkling  fail  simultaneously 


a  c: 
yi 


q  -  - 

iy  2 ,  . 


ca 


— >  c  - 


c  c 
yt 


accip  /p  ) 

c  • 


2A/3 


”  2  .  . 


y  t  3/2t. 

•>  p  /p  -  (  —  ) 

c  t 

o 


Therefore,  the  weight  function  becomes 
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w  -  Ip^aht  +  pabc  (  or  W  -  +  p  ) 


-  2p  ab 


a  c 
yt 


y  t  3/2A 


+  p  (  -  )  ahc 


Letting  d^/bc  -  0,  we  find 


c 

Opt 


3  /  2  A 


P  <T 

1  y  f 


(  1  +  3 / 2A  ) 


1/2 

] 


C 

opt. 


a  c 
yi  opt 


(P  ) 
c  opt 


yf 


0.57C 

3  £  I 


3/tA 

)  -  constant 


From  the  failure  mode  map,  it  is  realized  that  there  i 
restriction  in  t/a  for  the  co-occurrence  of  face  yielding  and  i 
wrinkling  failures. 


c 

opt 


a 


-  <  _ 


(7  C 
yi  opt 


a 


4  ®3  f®4 


8. 


■(p  /P  ) 

C  1  op 


C 

y  s 


c 

y  I 


— > 


(b)  Face 


C  (.g  g  ) 

*  *3  £**  '  '  f 

<7  ^ - ; - <7  ^ - 

2  yi 


P  , 

»  y  I 


wrinkling-core  shearing  fail 


P  (6B-3)/2A 

yi 

(  -  ) 

Q 


simultaneously 


C  a  c(p  /p  )* 

4  y  t  c  % 


‘?c.  “ 


g  a 


s  a 

ace 
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— > 


t  l/B 


P  /P  ”  ( 

c  ( 


Co  c 
*  y » 


act(p  /p  ) 

e  • 


2A/  3 


"  2  .  . 


r  - 


ac 


2A/3B 

) 


The  weight  function  is 


W  -  2p  abr  +  p  abc 


Letting  dV/8c  -  0.  It  is  found 


p  (3*3B)a 


2p^(2A-3B>(ja^(g^+i/^f^) 


^  (2A+3>/3B  3B/(2A-eB*3> 

(  -  )  ) 

C  a 
*  yi 


(P  ) 

c 


opt 


p 

s 


<7a^^ 


Co  c 
«  y*  opt 


l/B 

) 


C 

opt 


o(p  /P  ) 

C  1 


2A  /  3 
opt 


C 

opt 


Similarly,  from  the  failure  mode  map,  there  is  a  rest 
in  p^/p  for  face  wrinkling -core  shearing  failure  mode, 
i .  e . 


(P  /P  ) 

C  • 


Opt 


Co  c 

*  y I  opt 


l/B 

)  s  p  /p 

c  « 


a 

yi 


a 


3/2A 

) 


Therefore , 
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(6B-2A-3)/2a 


2p  (2A-3B)C  "o  ^  (g  +1^  e  )  o 

t  4  yi  '^3  4  yi 

<5  i  - ; -  (  -  ) 

p  (3-3B)g  a  o 

i  f 


(c)  Face  yielding*core  shearing  fail  simultaneously 


c  ct 
yi 

*’!»  ”  2. 

«  is  g  ) 

— >  r  - 

a  c 
yi 

C  (p  /p  )®ff  c 
*  e  s  y « 


p  -  p 


e  t 


C  <7  C 

*  y  • 


1/E 

) 


The  weight  function  is 
V  -  2p  abt  +  p  abc 

r  c 

Letting  3V/ac  -  0.  We  can  find 


c 

Opt 


2p^B<7a*(gJ+^'^^^) 


P  (B-l)r  , 
1  yi 


C  (7 

*  y » 


1/B 

) 


B/(2E-1) 


C  ' 
opt 


<7« 


a 

yi 


c 

opt 


(P  ) 
c 


opt 


p 

1 


Co  c 

*  y ■  opt 


l/B 

) 


In  order  for  face  yielding  and  core  shearing  to  occur  at  the 
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same  time  there  is  a  limitation  in  according  to  the  failure 


mode  map.  That  is 


I/B 


y f  3/2A 


ip  /P)  -  ( 

c  %  opt 


-  )  >  p  /p  ~  {  —  ) 

C  <7  C  O 

*  yt  opt 


Therefore , 


2  2 


2p  BC  o  (g  g  )  o 

t  *  y»  '*3  yf  (68-3)/ 

q  >  - -  (  -  ) 


2A 


p  (B-1)<t  g  ‘ 

t  y 


(d)  Face  yielding  face  wrinkling  -core  shearing  fail 
simultaneously 


a  c  c 


<7,..  - 


ty  2 ,  . 


QCC(P  /P  ) 

c  t 


2A/  3 


"  2 


S  (g  g  ) 
®3  t®4 


C(P  /p  Vc  c 

4  c  t  y « 


<?  - 

cs 


The  three  equations  in  three  unknowns  give 


oaF  c 

’  ®e  yf  -3B/2A 


c  - 

opt 


C  a 

*  y» 


)  -  constant 


t  - 

opt 


-  constant 


a  c 
y  f  opt 
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a 

yi 

(p  )  “  p  (  -  )  -  constant 

C  opt  1 

a 

In  running  this  progran,  the  q  requirement  is  checked 
first.  The  optimal  design  is  evaluated  for  each  case  for  which  the 
condition  is  met.  Among  the  above  possible  failure  modes,  the  one 
with  the  minimum  weight  is  selected  as  the  overall  optimum  design. 


V.  Minimum  Weight  Design  for  Both  Stiffness  and  Strength 


In  practice,  the  minimum  weight  design  should  take  both 
stiffness  and  strength  requirements  into  account.  In  this 
program,  there  are  three  kinds  of  failure  mode  considered  for 
sandwich  plates. 

(a)  Stiffness  and  face  yielding  failure 


a  cc 
yi 


qa^ 

— >  c  -  - 

a  C 
yi 


qa 


A  - 


16D 


4S 


— > 


P  - 

c 


2 

*  e 


2  A 


4C  E  c  q 

(  • 


8E  tc 

t 


i;g 


) 


P 

t 
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Therefore,  the  weight  function  becomes 

2  ,  .  „  2  * 


fi^abqa  a  g ^(.g 


1/0 


W  -  2p  abt  + 
£ 


I 


yt 


AC  E  £7 

I  »  y  f 


-l/G 


yf 

Since  the  weight  of  the  sandwich  plate  is  a  real  number,  t'- ^ 
value  in  the  bracket  should  be  larger  than  zero.  From  this 
observation  it  is  concluded  that 


0  <  t  < 


(1-1/  ^)g  o  * 
£  *!>•£ 


There  is  a  constraint  for  face  yielding  failure. 


ip  /p  )  >  p  /p  -  ( 

C  I  opt  C  I 


y£  3/2A 


Within  the  range,  a  minimum  weight  can  be  found  by 
incrementally  increasing  t.  The  corresponding  face  and  core 
thicknesses  and  foam  density  are  the  optimal  design  values. 

(b)  Stiffness  and  face  wrinkling  failure 


QCC(p  /p  ) 

c  • 


2A/3 


’£w 


<7* 


— >  C  - 


ac 


ip  /P  ) 

c  • 


•2A/3 
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qa 


qa 


- +  -  g^ 

16D  ^  4S 


2  /-i  2x 

I  ) 


— >  C  - 


8AE  (g  g  ) 

t'®3  f** 


.  .  .2A/3  ^ 

(P  /P  )  + 

c  « 


4C  E  A 

t  ■ 


(P  /P  ) 

c  * 


-G 


The  weight  function  is 


W  -  2p  abc  +  p  abc 

i  e 


From  the  failure  mode  map,  it  is  known  that  there  is  a 
limitation  in  p  /p  for  occurrence  of  face  wrinkling. 

e  f 

a 


yt 


0  <  p  IP  <  ( 

e  t 


) 


3/2A 


Also,  there  is  another  constraint  for  face  wrinkling  failure. 

/-  /  .  \  /  /  vB-2A/3 


“S. 


The  value  which  corresponds  to  the  minimum  weight  can  be 
found  using  the  same  argument  as  that  for  the  constraints  of 
stiffness  and  the  face  yielding  failure  mode. 

(c)  Stiffness,  face  yielding  and  face  wrinkling 


9 


iy 


a  cc 
y£ 


a^ (g  +w  g  ) 
®3  r®* 


/  /  x2A/3 

QCC(P  /p  ) 

C  1 

2 


4  2 

qa  qa 

^ - g  +  -  g 

16D  4S  * 
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The  three  constraints  in  three  variables  can  be  solved  for  : 


(p  /P  ) 

C  1 


opt 


a 

yt 


a 


3/2A 

) 


constant 


2  2  2 
aoFCl-i/)  as  g  a 

f  '  ^  ®2  yj  -3G/2A 

c  -  -  +  — -  (  -  )  -  constant 

8AE  (g  +1/  p  )  AC  E  A  o 

f  3  I  4  4ft 


t 

opt 


o 

y  t 


c 

opt 


constant 


In  running  this  program,  the  three  failure  modes  considered 
above  will  be  compared  with  each  other  to  determine  the  overall 
mir.imu..i  weight  design. 

It  is  noted  that  the  stiffness-core  shearing  failure  mode  is 
not  included  in  this  program  because  it  ends  up  with  an 
impractically  large  core  thickness. 


2o 


.■  ;"oam-core;^,£^;g^  •  ■.  •  . 


Figure  4,  The  geonetry  of  a  circular  sandwich  plate 
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Appendix 


Simply  supported  circular  sandwich  plate  : 

3-ft/  7+3i^ 

i  t 

£■1  -  -  (r/a)^  +  £n(a/r)  ■  - 

1+w  4  ( 1-t-i/  ) 

f  '  f 

-  1  +  2£nCr/a) 


1  a 

g.  -  1/2  +  - in(r/a)  -  — - 

2  8r^ 


1  a 

g - £r.ir/a)  +  — 7 

2  6r" 


1  a‘ 

g. - -  - ir(r/a)  -  — - 

2  2  £r 


g, - in(r/a)  +  — - 

®  2  8r* 

S,  -  0 
-  1/2 
-  1/2 


Clamped  circular  sandwich  plate 
(1)  a  >  0.588r 

g^  -  (r/a)*  •  in(r/a) 

g^  -  1  +  2in{r/a) 


g.'g.-O 

g,  -  0 

g,’  g,-  1/2 
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(2)  a  <  0.586r 


3 


g  -  (r/a)  -  ln{r/a)  ■  - 

4 


-  1  +  2in(r/a) 

1  a' 

- in(r/a)  + 


'3  “5  ‘=’6 


8r 


-  0 

-  1/2 

Simply  supported  rectangular  sandwich  plate  ; 


(16)- 

6 


II 


(ni-l)/2  (n-l)/2 

(-1)  (-1) 

2 

mnn 


64 


^2  -  —  I  I 


mnn 


^3  ■  ■^6  — r  ^  ^ 


16  (-1)  (*1)  CS 
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g,-  —  II 


16  n 
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^,-  —  11 


bn 
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'e  3 
jr 


II 


,  .  .  (  n  -  1  )  /  2 
(-1)  a 


nn 


16 

^8  -  “7 


II 


(•1) 


(m-l  )  12 


mn 


2  2 

me  ^ 
where  n  -  -  +  n 


m,  n  -  1,  3,  5.  7, 


.23 
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PROGRAM  SVP 

************************************************************* 


*** 

***  THIS  PROGRAM  WAS  CREATED  BY  JONG -SHIN  HUANG  *** 

***  IN  AUGUST  1988  FOR  OPTIMUM  DESIGN  METHODS  OF  *** 

***  STRUCTURAL  SANDWICH  PANELS  *** 

★★★  ★★★ 


************************************************************** 

INTEGER  TY.CONFI.FAIL 

REAL  C5(10)  ,C6(10)  . EF , EC .GC  ,WS  ,  SS  .UDL.MCL.LCL.M. SD.D.YF, YS . C3  ,  A,-G8 
REAL  STEP. Z.SIGMAC.TAUC.SIGMAF.TAUF. AREA, DEFL,ES.ROS.CG,ROF,G,B,G9 
REAL  STIFF, WEIGHT. THETA. AL.RCP.NUF.G1.G2.AA. BE. P.C4,Q,PP.TL, GRIT 
REAL  OMEGA. S. CiaO)  ,C2(10)  ,RCS  .  DELTA.  STRl .  STR2  .  STR3  ,  G3 .  G4 ,  G5  ,  G6  ,  G7 
DOUBLE  PRECISION  TF.TC.ROC.COEFF.STREN.AFA 
OPEN (UNlT-7 . FILE- *  SWP . OUT  * . STATUS- ' NEW ' ) 

WRITE(*.l) 

1  FORMAT (T2. 'WHAT  KIND  OF  PROBLEM  DO  YOU  WANT  TO  RUN  ?',/.T10, 

$'  1  -  ANALYSIS  OF  SANDWICH  BEAMS' ./.TIO, 

S'  2  -  MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  IN  SANDWICH  BEAMS',/, 
STIO,'  3  -  MINIMUM  WEIGHT  DESIGN  FOR  STRENGTH  IN  SANDWICH  BL^uMS '  , 
$/,T10.'  4  -  MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  AND  STRENGTH  IN  SA 
SNDWICH  BEAMS', 

$/,T10,'  5  -  ANALYSIS  OF  SANDWICH  PLATES' ./.TIO, 

S'  6  -  MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  IN  SANDWICH  PLATES',/, 
STIO,'  7  -  MINIMUM  WEIGHT  DESIGN  FOR  STRENGTH  IN  SANDWICH  PLATES'. 
$/,T10,'  6  -  MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  AND  STRENGTH  IN  SA 
SNDWICH  PUTES') 

READ(*,*)TY 

IF  (TY.EQ.l)  WRITE(7,51) 

IF  (TY.EQ.2)  WRITE(7,52) 

IF  (TY.EQ.3)  WRITE(7,53) 

IF  (TY.EQ.4)  WRITE(7,54) 

IF  (TY.EQ.5)  WRITE(7,55) 

IF  {TY.EQ.6)  WRITE(7,56) 

IF  (TY.EQ.7)  WRITE(7.57) 

IF  (TY.EQ.S)  WRITE(7,58) 

51  FORMAT (//.TIO, 'ANALYSIS  OF  SANDWICH  BEAMS',/) 

52  FORMATC//, TIO, 'MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  IN  SANDWICH  BEA 
SMS' ./) 

53  FORMAT (//.TIO, 'MINIMUM  WEIGHT  DESIGN  FOR  STRENGTH  IN  SANDWICH  BEAM 
SS' ./) 

54  FORMAT (//.TIO, 'MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  AND  STRENGTH  IN 
$  SANDWICH  BEAMS',/) 

55  FORMAT (//.TIO, 'ANALYSIS  OF  SANDWICH  PLATES',/) 

56  FORMAT (//.TIO, 'MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  IN  SANDWICH  PLA 
$TES ' ,/) 

57  FORMAT (//.TIO, 'MINIMUM  WEIGHT  DESIGN  FOR  STRENGTH  IN  SANDWICH  PLAT 
$ES'./) 

58  FORMAT (//.TIO, 'MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  AND  STRENGTH  IN 
$  SANDWICH  PLATES',/) 


-*♦  INPUT  THE  MATERIAL  PROPERTIES  *** 

WRITE(*,*)' INPUT  THE  ELASTIC  MODULUS  OF  FACE  MATERIAL  (  ksi  )' 
READ(*.*)EF 

WRITE(*,*)' INPUT  THE  YIELD  STRENGTH  OF  FACE  MATERIAL  (  ksi  )' 
READ(*,*)YF 


WRITE(*.*)'INPl'T 
READ(*.*)ROF 
VRITE(*,*) 'INPUT 
READ(*.*)ES 
URITE(*,*)' INPUT 
READ(*.*)YS 
WRITEC*,*) 'INPUT 
READ(*.*)ROS 
WRITE(*,*)’ INPUT 
$0F  FOAM  CORE’ 
READ(*,*)C3 
WRITE (*,*) 'INPUT 
$RE' 

READ(*,*)A 
WRITE(*.*) 'INPUT 
$  FOAM  CORE' 
READ(*,*)CG 
WRITE(*.*) 'INPUT 
$' 

READ(*.*)G 
WRITE (*.*) ' INPUT 
$F  FOAM  CORE’ 
READ(*.*)C4 


THE  MASS  DENSITY  OF  FACE  MATERIAL  (  pcf  )' 

THE  ELASTIC  MODULUS  OF  SOLID  FOAM  (  ksi  )' 

THE  YIELD  STRENGTH  OF  SOLID  FOAM  (  ksi  )' 

THE  MASS  DENSITY'  OF  SOLID  FOAM  (  pcf  ) ' 

THE  PROPORTIONALITY  CONSTANT  FOR  ELASTIC  MODULUS 

THE  POWER  CONSTANT  FOR  ELASTIC  MODULUS  OF  FOAM  .CO 

THE  PROPORTIONALITY  CONSTANT  FOR  SHEAR  MODULUS  OF 

THE  POWER  CONSTANT  FOR  SHEAR  MODULUS  OF  FOAM  CORE 

THE  PROPORTIONALITY'  CONSTANT  FOR  SHEAR  STRENGTH  0 


WRITE (*,*) 'INPUT  THE  POWER  CONSTANT  FOR  SHEAR  STRENGTH  OF  FOAM  COR 
$E' 

READ(*.*)B 

WRITE(7.75)EF.YT,ROF,ES.YS,ROS.C3.A.CG.G.C4.B 
75  FORMAT (TIO, 'THE  ELASTIC  MODULUS  OF  FACE  MATERIAL  . T50 , F12 . 3 . T65  , 

$'(  ksi  )',//, TIO, 'THE  YIELD  STRENGTH  OF  FACE  MATERIAL  -',T50,F12.3 
$,T65,'(  ksi  )',//, TIO, 'THE  MASS  DENSITY  OF  FACE  MATERIAL  -'.TSO.Fl 
$2.3,T65,'(  pcf  )',//, TIO, 'THE  ELASTIC  MODULUS  OF  SOLID  FOAM  -' ,T50 
$,F12.3,T65, ' (  ksi  )',//, TIO, 'THE  YIELD  STRENGTH  OF  SOLID  FOAM  -' , 
$T50,F12.3,T65. ' (  ksi  )',//, TIO, 'THE  MASS  DENSITY  OF  SOLID  FOAM 
ST50,F12.3,T65. ' (  pcf  )',//, TIO, 'THE  PROPORTIONALITY  CONSTANT',/, 
STIO.'FOR  ELASTIC  MODULUS  OF  FOAM  CORE  -' ,T50 . F12 . 3 . //, TIO . 

S' THE  POWER  CONSTANT' ,/, TIO, 'FOR  ELASTIC  MODULUS  OF  FOAM  CORE  -'  , 
$T50,F12.3,//,T10, 'THE  PROPORTIONALITY  CONSTANT' ./.TIO, 

S'FOR  SHEAR  MODULUS  OF  FOAM  CORE  -' ,T50 . F12 . 3 ,//,T10 . 

S’ THE  POWER  CONSTANT' ,/, TIO, 'FOR  SEAR  MODULUS  OF  FOAM  CORE  -' , 

ST50 , F12 . 3 , // , TIO . ' THE  PROPORTIONALITY  CONSTANT ' . / , TIO , 

S'FOR  SHEAR  STRENGTH  OF  FOAM  CORE  .T50 , F12 . 3 .//. TIO , 

S' THE  POWER  CONSTANT' ,/, TIO, 'FOR  SHEAR  STRENGTH  OF  FOAM  CORE  -' , 
ST50.F12.3./) 

AFA-O .57*C3**(2./3.) *ZF** ( 1 . /3 . ) *ES** ( 2 . /3 . ) 

GO  TO  (100, 100, 100, 100, 300, 300, 300, 300), TY 


■*"*-**it-*Tt-****^HHt’*  *•*-*■**  *  hlfk  •A-***'*-*-*******'*-'*-*****  <r*****^KnHr 
♦■y******  **  ■*■*♦<:**♦*' *’****-*'*T»***************'*St-  HitlrliHi 

**  A.  ANALYSIS  AND  DESIGN  OF  SANDWICH  BEAMS  *** 

■»■♦★****♦♦★★♦★★♦★*♦★★*★***★***********★********★★***★*★★★* 
★  ★★♦★★^^^^►♦★★^^^►*i^******ilrt^******Tlr**********'*T*-*-**********-* 


•****★*■*•*♦**★*******-*■★*********★****■*-******■*-*-*■**■*■  ******★*★★*•*•* 

**  INPUT  THE  CONFIGURATION  AND  LOADING  GEOMETRY  *** 

100  DATA  C1(1).C1(2).C1(3).C1(4).C1(5).C1(6)/4.1.8.8.2.12/ 

DATA  C2(1).C2(2).C2(3).C2(4),C2(5).C2(6)/2,1.2.2,1,2/ 

DATA  C5(1),C5(2) ,C5(3) ,C5(4),C5(5) ,C5(6)/48.3.192.76.fc.8,384/ 
DATA  r6(l).C6(2).C6(3),C6(4).C6(5),C6(6)/4.1.4.8.2.8/ 
WRITE(*.101) 


101  FORMAT(T2.’  WHAT  IS  THE  CONFIGURATION  AND  LOADING  GEOMETRY 
:  $T15,'  1  -  SIMPLY- SUPPORTED  UNDER  CONCENTRATED  LOAD',/. 

ST15,'  2  -  CANTILEVER  UNDER  CONCENTRATED  LOAD',/, 

$T15,'  3  -  FIXED  ENDS  UNDER  CONCENTRATED  LOAD',/, 

•  $T15.'  U  -  SIMPLY- SUPPORTED  UNDER  UNIFORM  DISTRIBUTED  LOAD’ ,/, 
$T15,'  5  -  CANTILEVER  UNDER  UNIFORM  DISTRIBUTED  LOAD',/, 

$T15.'  6  -  FIXED  ENDS  UNDER  UNIFORM  DISTRIBUTED  LOAD’) 
READ(*,*)CONFI 
IF  (CONFI.EQ.l)  WRITE(7,111) 

IF  (CONFI,EQ.2)  WRITE(7.112) 

IF  (CONFI.EQ.3)  WRITE (7, 113) 

IF  (C0NFI.EQ.4)  WRITE(7,114) 

IF  (CONFI.EQ.5)  WRITE(7.115) 

IF  (CONFL.EQ.6)  WRITE(7,116) 

111  FORMAKTIO, 'SIMPLY-SUPPORTED  UNDER  CONCENTRATED  LOAD',/) 

112  FORMATCTIO, 'CANTILEVER  UNDER  CONCENTRATED  LOAD',/) 

113  FORMAT '710. 'FIXED  ENDS  I'NDER  CONCENTRATED  LOAD'./) 

114  F0RMAT(T10, 'SIMPLY-SUPPORTED  UNDER  UNIFORM  DISTRIBUTED  LOAD',/) 

115  FORMAT (TIO, 'CANTILEVER  UNDER  UNIFORM  DISTRIBUTED  LOAD',/) 

116  FORMAT (TIO, 'FIXED  ENDS  UNDER  UNIFORM  DISTRIBUTED  LOAD',/) 
WRITE(*.*) 'INPUT  THE  WIDTH  OF  THE  SANDWICH  BEAM  (in.  )' 
READC*.*)WS 

WRITE (*,*) 'INPUT  THE  SPAN  OF  THE  SANDWICH  BEAM  (in.  )' 

READ(*,*)SS 

WRITE(7.117)WS,SS 

117  FORMAT(T10,  "THE  WIDTH  OF  THE  SANDWICH  BEAM  -' .T50 , FIO . 3 .165 , ' (  in. 
$  )',//, TIO, 'THE  SPAN  OF  THE  SANDWICH  BEAM  -' .T50 , FIO . 3 . T65 . ' (  in. 
$  )',/) 

IF  (TY.EQ.2)  GO  TO  200 
IF  (TY.EQ.3)  GO  TO  250 
IF  (TY.EQ.4)  GO  TO  270 

*********************************************** 

*********************************************** 

***  1,  analysis  of  SANDWICH  BEAMS  *** 

*********************************************** 
*********************************************** 

WRITE(*.*) 'INPUT  THE  MASS  DENSITY  OF  FOAM  CORE  (  pcf  )’ 
READ(*.*)ROC 

WRITE(*.*) 'INPUI  'THE  THICKNESS  OF  FACE  MATERIAL  (  in.  )’ 
READ(*.*)TF 

WRITE(*,*) 'INPUT  THE  THICKNESS  OF  FOAM  CORE  (in.  )' 

READ(*.*)TC 
WRITE(7,125)ROC.TF,TC 
EC-C3* (ROC/ROS ) **A*ES 
GC-CG* ( ROC/ROS ) **G*ES 
IF  (C0NFI.GE.4)  THEN 

WRITE(*,*) 'INPUT  THE  MAGNITUDE  OF  UNIFORM  LOAD  (  kips/in  )' 
READ(*.*)UDL 
WRITE(7,126)UDL 
ELSE  IF  (CONFI.lt. 4)  THEN 

WRITE(*,*) 'INPUT  THE  MAGNITUDE  OF  THE  CONCENTRATED  LOAD  (  kip 

$s  )' 

READ(*,*)MCL 

WRITE(*,*)' INPUT  THE  LOCATION  OF  THE  CONCENTRATED  LOAD  (  in. 

$)' 

WRITE (*,*)'***  THE  LOCATION  IS  MEASUF.ED  FROM  THE  FREE  END  FOR 
$CANTILEVER  BEAM  *★*' 

READ(*.*)LCL 

WRITE(7.127)MCL.LCL 


AL-(SS-LCL)/SS 

■  IF  (AL.lt. 0.01)  WRITE(*.*) '**  CHECK  THE  LOCATION  OF  LOAD  **' 
END  IF 

125  FORMAT (TIO, 'THE  MASS  DENSITY  OF  TORE  MATERIAL  .T50 ,F12 . 3 .T65 . 

•  $'(  pcf  )' .//.TIO.'THE  THICKNESS  OF  FACE  MATERIAL  .T50 . F12 . 3 . 

$T65.'(  in.  )',//, TIO. 'THE  THICKNESS  OF  FOAM  CORE  -' .T50 , F12 . 3 . 
$T65.'(  in.  )',/) 

126  FORMAT (TIO, 'THE  MAGNITUDE  OF  UNIFORM  LOAD  -' ,T50 . F12 . 3 . 

$T65.'(  kips/in.  )' ./) 

127  rr-MATCTlO.'THE  MAGNITUDE  OF  CONCENTRATED  LOAD  -' .T50 , F12 . 3 . 
$T65.'(  kips  )',//, 

$T10. 'THE  LOCATION  OF  CONCENTRATED  LOAD  -' .T50,F12.3, 

$T65.'(  in.  )',/) 


*★*  CALCULATE  THE  CRITICAL  BENDING  AND  SHEAR  STRESSES  *** 

IF  (CONFI.EQ.l)  THEN 

IF  (AL.lt. 0.5)  AL-I-AL 
C1(1)-1./(AL*(1-AL)) 

C2(1)*-1./AL 

C5(l)-27./(AL*(l-AL)*(2-AL)*(3*AL*(2-AL))**0.5) 

C6(1)-1./(AL*(1-AL)) 

ELSE  IF  (CONFI.EQ.2)  THEN 
C1(2)-1./AL 
C2(2)-l. 

C5 ( 2 ) -6 . / ( 3*AL*AL- AL**3 ) 

C6(2)-1./AL 

ELSE  IF  (CONFI.EQ.3)  THEN 
IF  (AL.CT.0.5)  AL-l-AL 
C1(3)-1./(AL*(1-AL)**2) 

C2 ( 3 ) -1 . / ( ( 1 - AL) **2* ( 1+2*AL) ) 

C5 ( 3 ) -3* ( 3  - 2*AL) **2/ ( 2*AL**2* ( 1 - AL) **3 ) 

C6 ( 3 ) -1 . / ( AL* ( 1 - AL) **2* ( 1+2*AL) ) 

END  IF 
SD-TF+TC 

D-EF*WS*TF**3/6 .  +EF*V’S*TF*SD**2/2  .  +EC*WS*TC**3/12 . 

IF  (C0NFI.CE.4)  THEN 

M-UDL’^SS**2/C1  (  CONFI ) 

Q-UDL*SS/C2(C0NFI) 

ELSE  IF  (CONFI.lt. 4)  THEN 
M-MCL*SS/C1( CONFI) 

Q-MCL/C2( CONFI) 

END  IF 
WRITE(*,130) 

WRITE(7.130) 

130  FORMAT (/,T5, 'THE  CRITICAL  BENDING  AND  SHEAR  STRESSES  TIO, 

$'2  (  in  )',T25,'SICMAC  (  ksi  )',T40,'TAUC  (  ksi  )',/) 
STEP-TC/40. 

Z-0.0  • 

DO  140  1-1,21 
SIGMAC-M*Z*EC/D 

TAUC-(EF*TF*SD/2 . +EC*(TC**2/4 . - Z*Z) /2 . )*Q/D 
WRITE(* , 135 ) Z , SICMAC , TAUC 
WRITE(7 , 135 ) Z . SIGMAC , TAUC 
135  FORMAT(T5,F7.3,T20,F12.6,T35,F12.6) 

Z-Z+STEP 
140  CONTINUE 
WRITE(*,145) 

WRITE(7,145) 


145  FORMAT  (/.  no, 'Z  '  in  )  *  ,T25 ,  '  SIGMAF  (  ksi  )'.T40.'TAUF  (  ksi  )',/) 
STFP-TF/20 . 

Z-TC/2 . 

DO  150  1-1,21 
SIGMAF-M*Z*EF/D 

TAUF-<5*EF*  (  (  (TC/2 .  +TF)  **2  -  Z*Z)  /2  . ) /D 
VT.ITE  (*,  135 )  Z .  SIGMAF .  TAUF 
WRITE ( 7 , 1 3 5 ) Z , S IGMAF , TAUF 
Z-Z+STEP 
150  CONTINUE 


**♦  CALCULATE  THE  MAXIMUM  DEFLECTION  *** 

AREA-WS*SD**2/TC 
IF  (C0NFI.GE.4)  THEN 

DEFL-UDL*SS**4/  (C5  ( CONFl )  *D )  +UDL*S $*’*•2/  ( AREA*GC*C6  ( CONFI ) ) 
ELSE  IF  (CONFI.lt. 4)  THEN 

DEFb-MCL*SS**3/ ( C5 ( CONFI ) *D ) +MCL*S  S / ( C6 ( CONFI ) *AREA*GC) 

END  IF 

WRITE(*,160)DEFL 
WRITE (7.1 60 )DEFL 

160  FORMAT (/.TIO, 'THE  MAXIMUM  DEFLECTION  ,T35 , ri2 . 6 ,T50 , ' (  in.  )') 

****************************************************** 

**★  JUDGE  THE  STRENGTH  FAILURE  MODE  *** 

PFT-Cl ( CONFI ) *YF*WS*TC*TF/SS 

PFW-AFA*C1 ( CONFI ) * (ROC/ROS )** ( 2 . *A/3 . ) *WS*TC*TF/SS 

COEFF- ( C3* (ROC/ROS ) ★*A*ES*SS/ ( 2 . *C1 ( CONFI ) *TF*EF ) ) **2 . 

COEFF-COEFF-*-  ( 1 .  /C2  ( CONFI ) )  **2 

PCS-C4* (ROC/ROS ) **B*YS*WS*TC* ( 1 . /COEFF) **0 . 5 

FAI1^3 

IF(PFy.LT.PFW)  THEN 

IF  (PFY.LT.PCS)  FAIl^l 
ELSE  IF  (PFW.LE.PFY)  THEN 
IF  (PFW.LT.PCS)  FAIL-2 
END  IF 

IF  (C0NFI.GE.4)  PP-UDL*SS 
IF  (C0NFI.lt. 4)  PP-MCL 
IF  (FAIL.EQ.l)  THEN 
WRITE(*,163) 

WRITE(7,163) 

IF  ((PFY-PP).GT.O)  THEN 
WRITE(*,161) 

WRITE(7,161) 

ELSE  IF  ((PFY-PP).LE.O)  THEN 
WRITE(*,162) 

WRITE (7. 162) 

END  IF 

ELSE  IF  (FAIL.EQ.2)  THEN 
WRITE (*,164) 

WRITE (7, 164) 

IF  ((PFW-PP) .GT.O)  THEN 
WRITE(*,161) 

WRITE(7,161) 

ELSE  IF  ((PFW-PP) .LE.O)  THEN 
WRITE(*,162) 

WRITE(7,162) 

END  IF 


ELSE  IF  (FAIL.EQ.3)  THEN 
.  WR1TE(*.165) 

WRITE(7,165) 

IF  ((PCS-PP).GT.O)  THEN 
WR1TE(*,161) 

WRITE(7,161) 

ELSE  IF  ((PCS-PP) .LE.O)  THEN 
WRITE(*,162) 

WRITE(7,162) 

END  IF 
END  IF 
CO  TO  1000 

161  FORMAKTIO, 'THE  SANDWICH  BEAM  UNDER  LOADING  IS  SAFE',/) 

162  FORMAT (TIO, 'THE  SANDWICH  BEAM  UNDER  LOADING  IS  NOT  SAFE',/) 

163  FORMAT (/.TIO, 'THE  FAILURE  MODE  IS  FACE  YIELDING',/) 

164  FORMAT(/,T10. 'THE  FAILURE  MODE  IS  FACE  WRINKLING',/) 

165  FORMAT (/, TIO, 'THE  FAILURE  MODE  IS  CORE  SHEAR  YIELDING',/) 


***  2.  MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  IN  SANDWICH  BEAMS  *** 

*★**★*♦**  *  *******■*-****■***********•*■•****★★★■**■»★********♦**  *  *  *  *  -kirirkirkie^trk-kifirieie 


***  INPUT  THE  REQUIRED  DESIGN  STIFFNESS 

200  WRITE(*,*)' INPUT  THE  REQUIRED  DESIGN  STIFFNESS  (  kips/in  )' 
READ(*,*) STIFF 

WRITE(7,201)STIFF 

201  FORMAT (TIO, ’THE  REQUIRED  DESIGN  STIFFNESS  -',T50, 
$F12.5,T65,'(  kips/in  )',/) 


■*  CALCULATE  THE  OPTIMAL  DESIGN  VALUES  *** 

C0EFF-4*(G*(4./(G-1) )**(!. /G)/(G-1))*((2+2*G)/(G-1))**(1-1./G) 
COEFF-COEFF* ( SS*C5 ( CONFI ) *EF/ ( CG*C6 ( CONFI ) *ES ))**(- 1 . /G ) 
COEFF-COEFF* ( STI FF/ ( C5 ( CONFI ) *WS*EF) ) ** ( 1 - 1 . /G ) 

TO  (  COEFF*ROP*SS**3/ROS  )  *★  ( G/  ( 3*G  - 1 )  ) 

TF-2* ( 1+G) *STIFF*SS**3/ ( (G - 1 ) *C5 ( CONFI ) *EF*WS*TC**2 ) 

COEFF-SS*C5 ( CONFI ) *STIFF*EF*TC*TF/ ( CG*C6 ( CONFI )*ES ) 

ROC- ( COEFF/ ( C5 ( CONFI ) *EF*WS*TF*TC**2 - 2*STIFF*SS**3 ) ) *★ (1 . /c ) *ROS 
WEIGHT-( 2 . *ROr*WS*TF*SS+ROC*WS*TC*SS )/1728000 . 

SD-TF+TC 

COEFF- ( CG* (ROC/ROS ) **g*ES*TC* ( 1+3* ( SD AF)**2 ) / ( 2*EF*TF) ) **0 . 5 
THETA-COEFF+SSAC 

WRITE (* , 215)TC , TF , ROC .WEIGHT . THETA 
WRITE(7 . 215 )TC . TF . ROC . WEIGHT , THETA 
215  FORMAT(//.T10.'THE  OPTIMAL  DESIGN  VALUES  :'.// 

$.T15.'THE  OPTIMAL  CORE  THICKNESS  -' ,T55.F12.4.T70. 'in. ' .//. 
$T15.'THE  OPTIMAL  FACE  THICKNESS  -' .T55 .F12.4.T70, ' in. ' .//.T15 , 
$'THE  OPTIMAL  MASS  DENSITY  OF  FOAM  -' ,T55,F12.4.T70. 'pcf ' ,//,T15, 
$'THE  MINIMUM  WEIGHT  OF  SANDWICH  BEAM  -' .T55.F12.4,T70. 'kips' ,//, 
§T15, 'THETA  OF  THE  SHEAR  LAG  CRITERION  -' .T55,F12.4) 

GO  TO  1000 


3.  MINIMUM  WEIGHT  DESIGN  FOR  STRENGTH  IN  SANDWICH  BEAMS 


*★*★^MHh*^^♦**•(M^*^H^**^^★★★*★★**^*-*****♦♦***★★★**★★★★★★★*****★*****♦****★★**★* 

*i.-ik*****^***Tlrtf»******************-ir*****’******T*-*'********************^**-*** 


***.  INPUT  THE  REQUIRED  DESIGN  STRENGTH  *** 

250  WRITE(*,*)' INPUT  THE  REQUIRED  DESIGN  STRENGTH  PER  UNIT  WIDTH  AND  L 
$ENCTH  (  p*i  )' 

READ(*,*)STREN 
WRITE(7,251)STREN 
STREN-STREN/1000 . 

251  FORMAT(T10.*THE  REQUIRED  DESIGN  STRENGTH  './.TIO, 

$'PER  UNIT  WIDTH  AND  LENGTH  ,T50.F12.4.T65. *psi' ,/> 


CALCULATE  THE  OPTIMAL  DESIGN  VALUES 


**★★**★  FY/rW/CS  FAILURE  DESIGN  ****** 

TC-STREN*SS*(YF/AFA)**( -1 . 5*B/A)/(YS*C4*C2(CONFI ) ) 

TF-STREN*S  S**2 . / ( Cl ( CONFI ) *YF*TC ) 

R0C-(YF/AFA)**(1 . 5/A)*R0S 

WEIGHT-(2 .*ROF*WS*TF*SS+ROC*WS*TC*SS)/1728000 . 

FAIL-4 

COEFF-(YF/AFA)**( (6*B-3)/(2.*A))*2*ROF*yS**2*C4**2*C2 (CONFI )**2 
STR1-C0EFF/(C1 ( CONFI ) *R0S*YF) 

COEFF- ( YF/AFA) ** ( ( 6*B - 2*A -  3 ) / ( 2 . *A ) ) * ( C2 ( CONFI ) *  C4*YS ) **2*R0F*2 
STR2-COEFr* ( 2*A - 3*B ) / ( ROS  * ( 3 . -  3 . *B ) *C1 ( CONFI ) *AFA) 

COEFF- (YF/AFA) **(( 6*B- 3 )/( 2 . *A) )*( C2 ( CONFI ) *C4*YS ) **2*B*2*ROF 
STR3-COEFF/ (ROS* ( B • 1 ) *C1 ( CONFI ) *YF) 

******  FW/FY  FAILURE  DESIGN  ****** 

IF  (STREN.LE.STRl)  THEN 

COEFF-2*ROr*STREN*SS**2*AFA** ( 1 . 5/A) /(Cl ( CONFI ) *R0S ) 

TCI- ( COEFF/ (YF** ( 1+ ( 1 . 5/A) ) ) ) **0 . 5 
TF1-STREN*SS**2/(C1 ( CONFI )*YF*TC1) 

ROCl- (YF/AFA) ** ( 1 . 5/A) *ROS 
W1-(2*TF1*WS*SS*ROF+TC1*WS*SS*ROC1)/1728000, 

IF  (Wl.LT, WEIGHT)  THEN 
TF-TFl 
TC-TCl 
ROC-ROCl 
WEIGHT-Vl 
FAIL-1 
END  IF 
END  IF 

*****  FW/CS  FAILURE  DESIGN  ****** 

IF  (STREN.LE.STR2)  THEN 

COEFF-<STREN*SS/(C2 (CONFI )*C4*YS) )**( (2*A+3)/( 3 . *B) )*AFA 
COEFF-COEFT*Cl ( CONFI ) *ROS * ( 3  >  3*B ) / ( 2 . *ROF* ( 2*A - 3*B ) ) 

TCI- ( COEFF/ ( STREN*SS**2 ) ) ** ( 3*B/ ( 2 . *A • 6 . *B+3 . ) ) 

ROCl- ( STREN*SS/ ( C2 ( CONFI ) *C4*YS*TC1 ))**( 1 . /B ) *ROS 
TFl-STREN*SS**2/(Cl(CONn)*AFA*(ROCl/ROS)**(2*A/3 . )*TC1) 
Wl-(2*Tn*WS*SS*ROF+TCl*WS*SS*ROCl)/1728000. 

IF  (Wl.LT. WEIGHT)  THEN 
TF-TFl 
TC-TCl 
ROC-ROCl 
WEICHT-Wl 
FAII^2 
END  IF 


END  IF 

******  FY/CS  FAILURE  DESIGN  ****** 

IF  ( SIREN. GE.STR3)  THEN 

C0EFF-2*B*R0F* ( STREN*S  S ) ** ( ( B - 1 ) /B ) *SS/ ( C 1 ( CONFI ) *YF) 

COEFF- ( C2 ( CONFI ) *C4*YS ) ** ( 1 . /B ) *COEFF/ (ROS* ( B - 1 ) ) 
TCl-COEFF** ( B/ ( 2 . *B - 1 . ) ) 

TFI-STREN*SS**2/(C1(C0NFI)*YF*TC1) 

ROC 1- ( STREN*S  S/ ( C2 ( CONFI ) *C4*YS*TC1 ) ) ** ( 1 . /B ) *ROS 
W1-(2*TF1*WS*SS*ROF+TC1*WS*SS*ROC1)/1728000. 

IF  (Wl.LT. WEIGHT)  THEN 
TF-TFl 
TC-TCl 
ROC-ROCl 
WEIGHT-Wl 
FAILp-3 
END  IF 
END  IF 

WRITE (*, 2 60 )TC,TF, ROC, WEIGHT 
WRITE ( 7 , 260 ) TC , TF , ROC , WEIGHT 
IF  (FAIL.EQ.l)  WRITE(7,261) 

IF  (FAIL.EQ.2)  WRITE(7,262) 

IF  (FAIL.EQ.3)  WRITE(7,263) 

IF  (FAIL.EQ.4)  WRITE(7,264) 

GO  TO  1000 

260  FORMAT(/,T10, 'THE  OPTIMAL  DESIGN  VALUES  :'.//,T15, 

$'THE  OPTIMAL  CORE  THICKNESS  ,T55 ,F12 .4,T70, ' in. ' ,//,T15 , 

$’THE  OPTIMAL  FACE  THICKNESS  .T55 , F12 .4,T70, ' in. ' ,//,T15 , 

$'THE  OPTIMAL  MASS  DENSITY  OF  FOAM  ,T55,F12.4,T70, 'pcf ' ,//,T15, 
$'THE  MINIMUM  WEIGHT  OF  SANDWICH  BEAM  ,T55,F12.4,T70, 'kips' ,/) 

261  FORMAT (T15,'FY/FW  FAILURE  DESIGN',/) 

262  FORMAT(T15, 'FW/CS  FAILURE  DESIGN',/) 

263  FORMAT (T1 5, 'FY/CS  FAILURE  DESIGN',/) 

264  FORMAT (T1 5, 'FY/FW/CS  FAILURE  DESIGN',/) 

******************************************************** 

******************************************************** 

***  4.  MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  *** 

***  and  STRENGTH  IN  SANDWICH  BEAMS  *** 

******************************************************** 
******************************************************** 

********************************************* 

***  INPUT  THE  DESIGN  PARAMETERS  *** 

********************************************* 

270  WR1TE(*,*)' INPUT  THE  REQUIRED  DESIGN  LOAD  (  kips  )' 

READ(*,*)P 

WRITE(*,*)' INPUT  THE  REQUIRED  DESIGN  DEFLECTION  (  in  ) ' 

READ (*,*) DELTA 
WRITE(7,271)P, DELTA 

271  FORMAT(T10, 'THE  REQUIRED  DESIGN  LOAD  -' ,T50,F12 .4,T65, 'kips' ,//, 
$T10,'THE  REQUIRED  DESIGN  DEFLECTION  -' ,T50,F12 .4,T65 , ' in. ' ,/) 

***************************************************** 

***  CALCULATE  THE  OPTIMAL  DESIGN  VALUES  *** 

***************************************************** 

***  S+FY  FAILURE  DESIGN  *** 

RCS-(YF/AFA)**(1 . 5/A)*ROS 

CRIT-DELTA*C5(CONFI)*EF*P/(2*C1<CONFI)**2*YF**2*WS*SS*1000 . ) 
TFl-CRIT 


WEIGHT-1000000 . 

•  FAll-1 

WRITE(*.*)'THE  PROGRAM  IS  RUNNING.  PLEASE  WAIT!' 

DO  273  1-1,999 

TC1-P*SS/ ( Cl( rONFI) *YF*WS*TF1 ) 

COEFF-C5 (CONFI ) *EF*P*SS*TF1*TC1/ ( CG*C6 (CONFI ) *£5 ) 
C0EFF-C0EFF/(C5(C0NFI)*DELTA*WS*TF1*TC1**2*EF-2*P*SS**3) 
ROCl-COEFP** ( 1 . /G ) *ROS 

W1-(2*ROF*WS*SS*TF1+ROC1*WS*SS*TC1)/1728000. 

IF  (ROCl.GE.RCS)  THEN 

IF  (Wl.LE. WEIGHT)  THEN 
TC-TCl 
TF-TFl 
ROC-ROCl 
WEIGHT-Wl 
END  IF 
END  IF 
TFl-TFl+CRlT 

273  CONTINUE 

***  S+FV  FAILURE  DESIGN  *** 

CRIT-RCS/1000. 

ROCI-CRIT 
DO  274  1-1.999 

C0EFF-2*SS**2*AFA*C1 ( CONFI ) * ( ROCl/ROS ) ** ( 2*A/3 . ) /C5 ( CONFI ) 
TC 1-COEFF/ ( DELTA*EF) +P*SS/ ( CG*C  6 ( CONFI ) *ES*WS*DELTA* 

$  (ROC1/ROS)**G) 

TF1-P*SS* (ROCl/ROS ) ** ( - 2*A/3 . ) / ( AFA*C1 ( CONFI ) *WS*TC1 ) 
W1-(2*ROF*WS*SS*TF1+ROC1*WS*SS*TC1)/1728000. 

TL-C2 ( CONFI ) *C4*ys* (ROCl/ROS ) ** ( B • 2*A/3 . ) / CAFA*C1 (CONFI ) ) 
IF  ((TF1/SS).LE.TL)  THEN 

IF  (Wl.LE, WEIGHT)  THEN 
TC-TCl 
TF-TFl 
ROC-ROCl 
WEIGHT-Wl 
FAII^2 
END  IF 
END  IF 

ROCl-ROCl+CRIT 

274  CONTINUE 

***  S+FY+FW  FAILURE  DESIGN  *** 

ROCl-(YF/AFA)**(l . 5/A)*ROS 

C0EFF-2*SS**2*C1 ( CONFI ) *YF/ ( C 5 ( CONFI ) *DELTA*EF) 
TCl-COEFF+P*SS*(Yr/AFA)**(-l . 5*G/A)/(CG*C6(CONFI)*ES*WS*DELTA) 
TF1-P*SS/ ( Cl ( CONFI ) ★Yr*WS*TCl ) 
W1-(2*ROF*WS*SS*TF1+ROC1*WS*SS*TC1)/1728000. 

IF  (Wl.LE. WEIGHT)  THEN 
TC-TCl 
TF-TFl 
ROC-ROCl 
WEIGHT-Wl 
FAIL-3 
END  IF 

WRITE (*. 2 80) TC , TF , ROC , WEIGHT 
WRITE(7,280)TC.TF, ROC. WEIGHT 
IF  (FAIL.EQ.I)  WRITE(7,281) 

IF  (FAIL.EQ.2)  WRITE(7,282) 


IF  (FAIL.EQ.3)  WRITE(7.283) 

GO,  TO  1000 

280  FORMAT (/.TIO, 'THE  OPTIMAL  DESIGN  VALUE  :',//,T15, 

$'THE  OPTIMAL  CORE  THICKNESS  .T55 . F12 .4 ,T70, ' in. ' T15 , 

.  $'THE  OPTIMAL  FACE  THICKNESS  .T55 .F12 .4 .T70, ' in. ' .//,T15 , 

$'THE  OPTIMAL  MASS  DENSITY  OF  FOAM  .T55 .F12 .4,T70 , 'pcf ' ,//.T15 . 

$*THE  MINIMUM  WEIGHT  OF  SANDWICH  BEAM  .T55.F12.4.T70, 'kips' ./) 

281  F0RMAT(T15.'ST1FFNESS/FY  FAILURE  DESIGN'./) 

282  FORMAT(T15, 'STIFFNESS/FV  FAILURE  DESIGN',/) 

283  F0RMAT(T15. 'STIFFNESS/FY/FW  FAILURE  DESIGN',/) 


it^eic 


B.  ANALYSIS  AND  DESIGN  OF  SANDWICH  PLATES 


***  INPUT  THE  CONFIGURATION  AND  LOADING  GEOMETRY  *** 

***  and  calculate  G1 - G9  COEFFICIENTS  **★ 

*****irk***1e***ir*ir******ieie***irieirir***ifleieit****ir*1c*irk***ir*-kirlHe*iriHt 

300  WRITE(*.301) 

301  FORMAT (T2. 'WHAT  IS  THE  CONFIGURATION  AND  LOADING  GEOMETRY  ?',/, 
$T15,'l  -  CIRCL'LAR  PLATE,  SIMPLY- SUPPORTED  UNDER  UNIFORM  LOAD.',/, 
$T15.'2  -  CIRCULAR  PLATE,  CLAMPED  UNDER  UNIFORM  LOAD.',/, 

$T15,'3  -  RECTANGULAR  PLATE,  SI  PLY -SUPPORTED  UNDER  UNIFORM  LOAD.') 

READ(*,*)C0NFI 
IF  (CONFI.EQ.l)  WRITE(7,311) 

IF  (CONFI.EQ.2)  WRITE<7,312) 

IF  (CONFI.EQ.3)  WP.ITE(7 , 313) 

311  FORMAT (TIO, 'CIRCULAR  PLATE,  SIMPLY- SUPPORTED  UNDER  UNIFORM  LOAD', 

$/) 

312  FOR.MAT  (TIO, 'CIRCULAR  PLATE.  CLAMPED  UNDER  UNIFORM  LOAD',/) 

313  FORMAT  (TIO, 'RECTANGULAR  PLATE.  SIMPLY -SUPPORTED  UNDER  UNIFORM  LOAD 

$'./) 

WRITE(*,*)' INPUT  THE  POISSON"  S  RATIO  OF  THE  FACE  MATERIAL' 

READ(*.*)NUF 

WRITE(7,323)NUF 

IF  (CONFI.LE.2)  THEN 

WRITE(*,*)' INPUT  THE  RADIUS  OF  THE  CIRCULAR  PLATE  (  in.  )' 

READ(*.*)RCP 

WRITE(7.321)RCP 

WRITE(*,*)' INPUT  THE  RADIUS  OF  THE  UNIFORM  LOAD  (  in.  )’ 

READ(*,*)AA 

WRITE(7.322)AA 

321  FORMAT (TIO, 'THE  RADIUS  OF  THE  CIRCULAR  PLATE  ,T50,F12.3,T65, 

S'(  in.  )',/) 

322  FORMAT (TIO, 'THE  RADIUS  OF  THE  UNIFORM  LOAD  ,T50,F12.3,T65, 

$'<  in.  )',/) 

323  FORMAT  (TIO, 'THE  POISSON"  S  RATIO  OF  THE  FACE  MATERIAL  .T50, 
$n2.3,/) 

IF  (CONFI.EQ.l)  THEN 

Cl-  (  3-f-NUF)  ★  (RCP/AA)  **2/  ( 1+NUF)  +L0G  (  AA/RCP ) 

$  -(7+3*NUF)/(4*( 1+NUF)) 

G3-0 . 5+0 . 5*LOG (RCP/AA) - AA**2/ ( 8+RCP++2 ) 

G4-0 , 5*LOG (RCP/AA)+AA**2/ ( 8*RCP**2 ) 

G5-G3 

G6-G4 


ELSE  IF  (C0NFI.EQ.2)  THEN 

G 1- ( RCP/AA ) **2 • LOG ( RCP/AA ) - 0 . 7 5 
IF  ((AA-0.588*RCP) ,GE.O)  THEN 
G3-0.5-AA**2/(4*RCP**2) 

G4-0. 

G5-;  . 

G6-G3 

ELSE  IF  ((AA-0.588*RCP) .LT.O)  THEN 

G3-0 . 5*LOG (RCP/AA )+AA**2/( 8*RCP**2 ) 
G4-G3 
G5-G3 
G6-G3 
END  IF 
END  IF 

G2-l+2*LOG (RCP/AA) 

G7-0. 

G8-0.5 

G9-0.5 


ELSE  IF  (CONFI.GE.3)  THEN 

WRITE (*,*) 'INPUT  THE  LENGTH  (  X  DIRECTION  )  OF  THE  RECT/UCG'JLA 
$R  PLATE  (  in.  )' 

READ(*.*)BB 

WRITE(7,324)BB 

WRITE (*,*) 'INPUT  THE  WIDTH  (  Y  DIRECTION  )  OF  THE  RECTANGULAR 
$PLATE  (  in.  )' 

READ(*,*)AA 
WRITE (7, 325 )AA 

324  FORMAT (TIO, 'THE  LENGTH  (  X  DIRECTION  )  OF',/.T15, 

$'THE  RECTANGULAR  PLATE  -' .T50 . F12 . 3 .T65 , ' (  in.  )',/) 

325  FORMAT (TIO, 'THE  WIDTH  (  Y  DIRECTION  )  OF',/,T15, 

§' THE  RECTANGULAR  PLATE  -' .T50,F12.3,T65, ' (  in.  )’,/) 

DO  330  1-0,11 
M-2*I+1 

DO  330  J-0.11 
N-2*J+1 


OMEGA- (M*AA/BB ) **2+N**2 

G1-G1+16**2*( • 1)**I* ( - l)**J/( 3 . 14159**6*M*N*0MEGA**2 ) 
G2-G2+64*(-l)**I*(-l)**J/(3.1415  9**4*M*N*0MEGA ) 

G3-C3+16* ( - 1 ) **I* ( - 1 ) **J*M*AA**2/ ( 3 . 1415  9**4*OMEGA**2*N* 
$  BB**2) 

G4-G4+16* ( - 1 ) **I* ( - 1 ) **J*N/ ( 3 . 14159**4*OMEGA**2*M) 
G7-G7+16*AA/(3 . 14159**4*BB*0MEGA**2) 

G8-G8+16* ( - 1 ) **J*AA/ (3.1415  9**3*N*OMEGA*BB ) 

G9-G9+16*( - l)**i/(3 . 14159**3*M*0MEGA) 

330  CONTINUE 
G5-G4 
G6-G3 
END  IF 

IF  (Ty.EQ.6)  GO  TO  400 
IF  (TY.EQ.7)  GO  TO  500 
IF  (TY.EQ.8)  GO  TO  600 


1.  ANALYSIS  OF  SANDWICH  PLATES 


WRITE(*,*)' INPUT  THE  MAGNITUDE  OF  THE  UNIFORM  LOAD  (  ksi  ) ' 
READ(*,*)UDL 

WRITE(*.*)' INPUT  THE  THICKNESS  OF  FACE  MATERIAL  (  in.  )' 


Ri>D(*,*)TF 

miTEC*,*)' INPUT  the  thickness  OF  FOAM  CORE  (in,  )' 

READ(*,*)TC 

WRITE (♦.*) 'INPUT  THE  MASS  DENSITY  OF  FOAM  CORE  (  pcf  ) ' 
READ(*,*)ROC 

WRITE(7,341)UDL,TF,TC,ROC 

341  FORMAT (TIO, 'THE  MAGNITUDE  OF  THE  UNIFORM  LOAD  , T50 . F12 . 5 .T65 , 

$'(  ksl  )'.//. TIO. 'THE  THICKNESS  OF  FACE  MATERIAL  -' ,T50 , F12 . 3 ,T65 , 
$'(  In.  )',//, TIO, 'THE  THICKNESS  OF  FOAM  CORE  -' ,T50 . F12 . 3 .T65 . 

$'(  In.  )',//, TIO, 'THE  MASS  DENSITY  OF  FOAM  CORE  -' ,T50 . F12 . 3 .T65 , 
$'(  pcf  )',/) 


SD-TC+TF 

I>-EF*TF*SD**2/  (  2 .  *  ( 1 .  -  NUF**2 ) ) 

GC-CG* (ROC/ROS ) ★*G*ES 
S-SD*GC 

DEFL-UDL*AA**4*G1/ ( 16*D) +UDL*AA**2*G2 / ( 4*S ) 

SIGKAX-UDL*AA**2*  (G3+NTJF*G4  )  /  (  SD*TF) 

S1GMAY-UDL*AA**2* ( G5+NUF*G6 ) / ( SD*TF) 

TAUXY-UDL*AA**2* ( 1 - NUF) *07/ ( SD*TF) 

TAUZX-UDL*AA*G8/SD 

TAUYZ-UDL*AA*G9/SD 

WRITE ( * . 346 ) S I GMAX , S I GMAY , TAUXY , TAU2X , TAUYZ , DEFL 
WRITE  (  7 . 346  )  SIGMAX ,  S I  GMAY .  TAL^XY ,  TAUZX  ,  TAUYZ ,  DEFL 
346  FORMAT(T10. 'THE  MAXIMUM  STRESSES  IN  THE  FACE  MATERIAL 

STIS.'THE  MAXIMUM  STRESS  SIGMAX  (SIGMAR)  ,T50,F12 .4,T65 , ' ksi' , // , 

$T15.'THE  MAXIMUM  STRESS  SIGMAY  (SIGMAT)  ,T50.F12.4,T65, 'ksi' ,//, 

$T15.'THE  MAXIMUM  STRESS  TAUXY  (TAURT)  -' .T50.F12.4,T63, 'ksi' ,//, 
$T10,-THE  MAXIMUM  STRESSES  IN  THE  FOAM  CORE 
$T15,'THE  MAXIMUM  STRESS  TAUZX  (TAU2R)  -' ,T50,F12 .4,T65 , 'ksi ' 
$T15,'THE  MAXIMUM  STRESS  TAUYZ  (TAUTZ)  -' ,T50. F12 . 4,T65 , 'ksi ' ,// , 
$T10,'THE  MAXIMUM  DEFLECTION  -' ,T50,F12.4,T65. ' in. ' ) 


JUDGE  THE  STRENGTH  FAILURE  MODE  *** 

?FY-YF*TC*TF/ (AA**2* (G3+NUF*G4 ) ) 

PFW-AFA*TC*TF* (ROC/ROS )**( 2 *A/3 . )/(AA**2*(G3+NUF*G4) ) 
PCS-C4* (ROC/ROS ) **B*YS*TC/ ( AA*G8 ) 

FAI1^3 

IF(PFY.LT.PFW)  THEN 

IF  (PFY.LT.PCS)  FAIW 
ELSE  IF  (PFW.LE.PFY)  THEN 
IF  (PFW.LT.PCS)  FAII^2 
END  IF 

IF  (FAIL.EQ.l)  THEN 
WRITE(*,363) 

WRITE(7.363) 

IF  ((PFy-DDL).GT.O)  THEN 
WRITE(*,361) 

WRITE(7,361) 

ELSE  IF  ((PrY.UDL).LE.O)  THEN 
WRITE(*,362) 

WRITE(7.362) 

END  IF 

ELSE  IF  (FAIL.EQ.2)  THEN 
WRITE(*,364) 


WRITE (7, 364) 

IF  ((PFW-UDL) .GT.O)  THEN 
WRITE(*,361) 

WRITE<7,361) 

ELSE  IF  ( (PFW-UDL). LE.O)  THEN 
WRITE(*.362) 

WRITE(7,362) 

END  IF 

ELSE  IF  (FAIL.EQ.3)  THEN 
WRITE(*,365) 

WRITE(7.365) 

IF  ( (PCS -UDL). GT.O)  THEN 
WRITE(*,361) 

WRITE(7.36X) 

ELSE  IF  ((PCS -UDL). LE.O)  THEN 
WRITE(*,362) 

WRITE (7. 362) 

END  IF 
END  IF 
GO  TO  1000 

361  FORMAT(T10. 'THE  SANDWICH  PLATE  UNDER  LOADING  IS  SAFE',/) 

362  FORMAT(T10,'THE  SANDWICH  PLATE  UNDER  LOADING  IS  NOT  SAFE',/) 

363  FORMAT (/.TIO, 'THE  FAILURE  MODE  IS  FACE  YIELDING',/) 

364  FORMAT (/.TIO, 'THE  FAILURE  MODE  IS  FACE  WRINKLING',/) 

365  FORMAT (/, TIO, 'THE  FAILURE  MODE  IS  CORE  SHEAR  YIELDING',/) 


2.  MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  IN  SANDWICH  PLATE  *** 


INPUT  THE  REQUIRED  DESIGN  STIFFNESS 


400  WRITE(*,*)' INPUT  THE  REQUIRED  DESIGN  STIFFNESS  (  ksi/in  )' 
READ(*,*) STIFF 

WRITE (7, 401) STIFF 

401  FORMAT (TIO, 'THE  REQUIRED  DESIGN  STIFFNESS  .T50 ,F12 . 5 ,T65 , 
$'(  ksi/tn  )',/) 


Ifkliiilfk'kiHeieirkitieiritlftiii'tcliieli  k)t1i*'k-tie'k1rie1r1e**irit  *  *  ***  *  *  it-**** 

*  CALCULATE  THE  OPTIMAL  DESIGN  VALUES  *** 

ROC-( (G+1) **2* ( G- 1 ) **2*G2**3*ROS*EF* ( STIFF*AA) **2/ 

S ( 256*CG**3*G* ( 1 -NUF**2 )*G1*R0F*ES**3 ) ) ♦* ( 1 . / ( 3*C- 1 . ) ) *ROS 

TC-(4**(G+1 ) *CG*G**G* (G+1 ) ** (G • 1 ) * ( 1 - NUF**2 ) **G*G1**G*R0F**G*ES* 

S ( STI FF*AA) ** (G - 1 ) / ( ( C - 1 ) ** ( 2*G ) *G2*R0S**G*EF**G ) ) *★ ( 1 . / ( 3*G - 1 . ) ) 
$*AA/2 . 

TF-((G**2.1)**(G+l)*(l-NUF<rtr2)**(G-l)*Gl**(G-l)*G2**2*ROS**(2*G)* 

$Er**(l-G)*(STlFF*AA)**(G+l)/(2**(13*G+l)*CG'W2*G**(2*G)*ROF**(2*G) 

S*ES**2))**(1./(3*G-1.))*4*AA 

IF  (CONFI.lt. 3)  WEICHT-(2*TF*ROF+TC*ROC)*3.14*RCP**2/1728000. 

IF  (CONn.GE.3)  WEICHT-(2*TF*ROF+TC*ROC)*AA*BB/1728000. 

WRITE (*. 402 ) TF , TC , ROC . WEIGHT 

WRITE ( 7 , 402 ) TF , TC , ROC , WEIGHT 
.02  FORMAT(//,T10.'THE  OPTIMAL  DESIGN  VALDES  :',//, 

$T15,'THE  OPTIMAL  THICKNESS  OF  FACE  MATERIAL  ,T55,F12 .4, 
$T70,*in.' ,//,T15,'THE  OPTIMAL  THICKNESS  OF  CORE  MATERIAL  -'.T55, 
$F12.4,T70,'ln.',//,T15,'THE  OPTIMAL  MASS  DENSITY  OF  FOAM  -',T55. 


$F12.4,T70. 'pcf' ,//,T15. 'THE  MINIMUM  WEIGHT  OF  SANDWICH  PLATE  . 
$Tb5-.F12.4.T70,  'kips'  ./) 

GO  TO  1000 


*4^ 


3.  MINIMUM  WEIGHT  DESIGN  FOR  STRENGTH  IN  SANDWICH  BEAMS  *** 


***  INPUT  THE  REQUIRED  DESIGN  STRENGTH 

500  WRITE(*.*)' INPUT  THE  REQUIRED  UNIFORM  LOADING  (  psi  )' 
READ(*,*)STREN 
WRITE(7,551)STREN 
STREN-STREN/1000. 

551  FORMAT (TIO, 'THE  REQUIRED  UNIFORM  LOADING  -' .T50,F12 .4 ,T65 , 'psi' ,/) 


CALCULATE  THE  OPTIMAL  DESIGN  VALUES 


*-** 


**-«r***  FY/FW/CS  FAILURE  DESIGN  *•**♦★* 

TC-STREN*AA*G8*(YF/AFA)*-*( -1 . 5*B/A)/(YS*C4) 

TF-STREN*AA**2 . * (G3+NUF+G4 ) / ( YF*TC ) 

R0C-(YF/AFA)**(1 . S/A)*ROS 

IF  (CONFI.lt. 3)  WEIGHT-(2*ROF*TF+ROC*TC)*3.14*RCP**2/1728000. 

IF  (C0NFI.GE.3)  WEIGHT- (2*ROF*TF+ROC*TC)*AA*BB/17 2 8000 . 

FAIb-4 

COEFF- (YF/AFA) ** ( ( 6*B • 3 ) / ( 2 . *A) ) ★2*ROF+YS**2*C4**2* ( G3+NUF*G4 ) 
STRl-C0EFF/(G8**2*R0S*yF) 

COEFF- (YF/AFA) ** ( ( 6*B - 2*A- 3 ) / ( 2 . *A) ) * (C4*YS ) **2*ROF*2* (G3+NUF*G4 ) 
STR2-COEFF*(2*A- :*B)/(ROS* ( 3 . -  3 . *B) *G8**2*aFA) 

COEFF- (YF/AFA)  **  (  (  6*B -  3 )  / ( 2  .  *A ) )  *  ( C4*YS  )  **2-^S*2*ROF*  ( G3+NUF*G4  ) 
STR3-COEFF/(ROS*(B-l)*G8**2*YF) 

******  FW/FY  FAILURE  DESIGN  ****** 

IF  (STREN.LE.STRl)  THEN 

COEFF-2*ROF*STREN*AA**2*AFA**(l . 5/A)*(G3+NUF*G4)/R0S 
TCl-( COEFF/ ( YF** ( 1+ ( 1 . 5/A) ) ) )**0 . 5 
TF1-STREN*AA**2* ( G3+NUF*G4 ) / (  YF*TCl ) 

R0C1-(YF/AFA)**(1 . 5/A)*R0S 

IF  (CONFI.lt. 3)  W1-(2*TF1*ROF+TC1*ROC1)*3.14*RCP**2/1728000. 
IF  (C0NFI.GE.3)  W1-(2*ROF*TF1+ROC1*TC1)*AA*BB/1728000 . 

IF  (Wl.LT. WEIGHT)  THEN 
TF-TFl 
TC-TCl 
ROC-ROCl 
WEIGHTS 
FAII^l 
END  IF 
END  IF  ■ 

******  FW/CS  FAILURE  DESIGN  ****** 

IF  (STREN.LE.STR2)  THEN 

C0EFF-( STREN*AA*G8/ ( C4*YS ))**(( 2*A+3 ) / ( 3 . *B ) ) *AFA 
C0EFF-C0EFF*R0S* ( 3 - 3*B) / ( 2 . *ROF* ( 2*A- 3*B) ) 

TCI- ( COEFF/ ( STREN*AA**2* ( G3+NUF*G4 ) ) ) ** ( 3*B/ ( 2 . *A - 6 . *B+3 . ) ) 
R0C1-(STREN*AA*G8/(C4*YS*TC1) )**(!. /B)*ROS 
TFl-STREN*AA**2*(C3+NUF*G4)/(AFA*(R0Cl/R0S)**(2*A/3 . )*TC1) 

IF  (CONFI.lt. 3)  W1-(2*TF1*ROF+TC1*ROC1)*3.14*RCP**2/1728000. 
IF  (C0NFI.GE.3)  W1-(2*TF1*ROF+TC1*ROC1)*AA*BB/1728000 . 

IF  (Wl.LT. WEIGHT)  THEN 


ROC-ROCl 
WEIGHT-Vl 
FAI1^2 
END  IF 
END  IF 

******  fY/CS  faiujre  design  ****** 

IF  ( SIREN. GE.STR3)  THEN 

C0EFF-2*B*R0F* ( STREN*AA) ** ( ( B • 1 ) /B ) *AA* (G3+NUr*G4 ) /YF 
TCl-((C4*YS/G8)**(l./B)*COEFF/(ROS*(B-l)))**(B/(2,*B-l.)) 
TF1-STREN*AA**2* ( G3+NUF*C4 ) / ( Yr*TCl ) 
R0C1-(STREN*AA*G8/(C4*YS*TC1) )**(! ./B)*ROS 
IF  (CONTI. LT. 3)  W1-(2*TF1*R0F+TC1*R0C1)*3 . 14*RCP**2/1728000 
IF  (CONFI.GE.3)  Wl-(2*TFl*R0F+TCl*ROCl)*AA*BB/1728000 . 

IF  (Wl.LT. WEIGHT)  THEN 
TF-TFl 
TC-TCl 
ROC-ROCl 
WEIGHT-Vl 
FAII^3 
END  IF 
END  IF 

WRITE(* . 560)TC , TF. ROC . WEIGHT 
WRITE ( 7 , 5 60 ) TC , TF . ROC , WEI GHT 
IF  (FAIL.EQ.l)  WRITE(7,561) 

IF  (FAIL.EQ.2)  WRITE(7,562) 

IF  (FAIL.EQ.3)  WRITE(7,563) 

IF  (FAIL.EQ.4)  WRITE(7,564) 

GO  TO  1000 

560  FORMAT(/,T10,'THE  OPTIMAL  DESIGN  VALUES  :'.//,T15, 

$'THE  OPTIMAL  CORE  THICKNESS  -',155 ,F12 .4,T70, ' in. ' ,//,T15 , 

S'THE  OPTIMAL  FACE  THICKNESS  ,T55,F12.4,T70, 'in. ' ,//,T15, 

S'THE  OPTIMAL  MASS  DENSITY  OF  FOAM  .T55 . F12 .4,T70, 'pcf ' . T15 , 
$'THE  MINIMUM  WEIGHT  OF  SANDWICH  PLATE  ,T55 .F12 .4,T70 , 'kips ’ ,/) 

561  FORMAT (TIS.'FY/FW  FAILURE  DESIGN',/) 

562  FORMAT (T15, 'FW/CS  FAILURE  DESIGN',/) 

563  FORMAT (T15, 'FY/CS  FAILURE  DESIGN',/) 

564  FORMAT (T15, 'FY/FW/CS  FAILURE  DESIGN',/) 


***  4.  MINIMUM  WEIGHT  DESIGN  FOR  STIFFNESS  *** 
***  AND  STRENGTH  IN  SANDWICH  PLATES  *** 
★  0  *********■************<*****************  *************** 

*********************************************************** 


********************************************* 

***  INPUT  THE  DESIGN  PARAMETERS  *** 

********************************************* 

600  WRITE(*,*)' INPUT  THE  REQUIRED  UNIFORM  LOADING  (  ksi  )' 

READ(*,*)P 

WRITE(*,*)' INPUT  THE  DESIGN  DEFLECTION  (  in  ) ' 

READ(*,*)DELTA 
WRITE(7,671)P, DELTA 

671  FORMAT (TIO, 'THE  REQUIRED  UNIFORM  LOADING  -' ,T50,F12.4,T65, 'kips' 
$//,T10,'THE  DESIGN  DEFLECTION  -' ,T50,F12.4,T65, 'in. ' ,/) 

*****  ******************  ******  *★*★★★*  *****  *■<  *  *  *  *  *■**★*■* 

**  CALCULATE  THE  OPTIMAL  DESIGN  VALUES 


»* *★★♦*** ■#r***^  ii***ic***irie*1r1rit***ir*ieirirk**ir**** 


***  S+nr  FAILURE  DESIGN  *★* 

.  RCS-(YF/AFA)**(1.5/A)*R0S 

CRIT-8*DELTA*EF*P*(C3+NUr*G4)**2/( (1-NUF**2)*G1*YF**2*1000 . ) 

TFl-CRIT 

WEIGHT-1000000. 

FAIb-1 

WRITE(*.'*-)*THE  PROGRAM  IS  RUNNING.  PLEASE  WAIT!' 

DO  673  1-1,999 

TC1-P*AA**2* (G3+NUF*G4 ) / ( YF*TF1 ) 

COEFF-(  (DELTA/P)  - AA*v-4»  ( 1 . NUF**2 )*Gl/( 8*EF*TF1*TC1**2 ))**(-.  1) 
R0C1-(C0EFF^AA**2*G2/(4*CG*TC1*ES) )**(! . /G)*ROS 
IF  (CONFI.lt. 3)  W1-(2*ROF*-TF1+ROC1*TC1)*3.14*RCP**2/1728000. 
IF  (C0NFI.GE.3)  W1-(2*ROF*TF1+ROC1*TC1)*AA*BB/1728000 . 

IF  (ROCl.GE.RCS)  THEN 

IF  (Ul.LE. WEIGHT)  THEN 
TC-TCl 
TF-TFl 
ROC-ROCl 
WEIGHT-Wl 
END  IF 
END  IF 

TFl-TFl+CRIT 

673  CON’TINUE 

***  S+FW  FAILURE  DESIGN  *** 

CRIT-RCS/1000 . 

ROCl-CRIT 
DO  674  1-1,999 

C0EFF-(R0C1/R0S )** ( 2*A/3 . ) *AA**2*AFA*G I* ( 1 -NUF**2 ) /8 . 
TCl-COEFF/ (DELTA*EF* (G3+NUF*G4 ) ) +P*AA**2*G2/ ( CG*4 . *ES*DELTA* 

$  (R0C1/R0S)**G) 

TF1-P*AA**2* (G3+NUF*G4 )* (ROCl/ROS )★*( - 2*A/3 . ) / ( ArA*TCl ) 

IF  (CONFI.lt. 3)  W1-(2*ROF*TF1+ROC1*TC1)*3.14*RCP**2/1728000. 
IF  (C0NFI.GE.3)  W1-(2*ROF*TF1+ROC1*TC1)*AA*BB/1728000 . 
Tl^C4*YS*(G3+NUF*G4)*(R0Cl/R0S)**(B-2*A/3 . )/(AFA*G8) 

IF  ((TF1/AA).LE.TL)  THEN 

IF  (VI. LE. WEIGHT)  THEN 
IC-TCl 
TF-TFl 
ROC-ROCl 
WEIGHT-Wl 
rAIL-2 
END  IF 
END  IF 

ROCl-ROCl+CRIT 

674  CONTINUE 

***  S+FY+FW  FAILURE  DESIGN  *** 

ROCl- (YF/AFA) ** ( 1 . 5/A ) ♦ROS 

C0EFF-AA**2*YF*G1*  ( 1  -  NUF<^*2  )  /  (  8*DELTA*EF*  (G3+NUF*G4  )  ) 
TC1-C0EFF+P*AA**2*G2*(YF/AFA)**( *1 . 5*G/A)/(4*CG*ES*DELTA) 
TF1-P*AA**2* (G3+NUF*G4 ) / ( YP*TC1 ) 

IF  (CONFI.lt. 3)  V1-(2*ROF*TF1+ROC1*TC1)*3.14*RCP**2/1728000. 

IF  (C0NFI.CE.3)  W1-(2*ROF*TF1+ROC1*TC1)*AA*BB/1728000. 

IF  (Wl.LE. WEIGHT)  THEN 
TC-TCl 
TF-TFl 
ROC-ROCl 


WEIGHT-Wl 
FAIL-3 
END  IF 

WRITE(* . 680)TC , TF , ROC . WEIGHT 
WRITE ( 7 . 6 80 ) TC . TF . ROC , WEIGHT 
IF  (FAIL.EQ.l)  WRITE(7,681) 

IF  (FAIL.EQ.2)  WRITE(7,682) 

IF  (FAIL.EQ.3)  WR1TE(7,683) 

680  FORMAT(/.T10/THE  OPTIMAL  DESIGN  VAUJE  :',//.Tl5. 

$'THE  OPTIMAL  CORE  THICKNESS  -* ,T55.F12.4,T70, 'in. ' ,//,Tl5, 

$'THE  OPTIMAL  FACE  THICKNESS  ,T55,F12.4.T70. 'in. ' .//,T15, 

$'THE  OPTIMAL  MASS  DENSITY  OF  FOAM  -' ,T55 ,F12 .4,T70, 'pcf ' ,//, T15 , 
$'THE  MINIMUM  WEIGHT  OF  SANDWICH  PLATE  -' .T55,F12.4.T70, 'kips' ,/) 

681  PORMAT(T15,'STIFFNESS/FY  FAILURE  DESIGN',/) 

682  FORMAT (T15.'STIFFNESS/FW  FAILURE  DESIGN',/) 

683  FORMAT(T15,'STIFFNESS/FY/FW  FAILURE  DESIGN',/) 


1000  STOP 
END 
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THE  MAXirr.lH  DEFLECT  I  Of-i  =  .035319  in.  ) 

THE  FA  I  M  IRE:  nOOE.  3E  F-.-.CE  ViEi  rCI  'F 

THE  SAI'lL'IfCH  beam  UM'-  -F!  IS  t>AFE 


ni  N I  Mill  1  tJrlGMr  Dt-- J I  .'I  J  FLU-;  STPiE-fUGl  H  ll'  ivir  JIH'J  ]  CH  E'EAMS 

THE  ELASTIC  MGr-ULliS  OF  FACE  MATE  RIAL  =  10X50.000  (  Lei  ) 

THE  YIELD  STRENGTH  OF  FACE  MATERIAL  =  lE.R7u  ^  Lei  ) 

THE  MASS  DENSITY  OF  FACE  nATE.F<IAL  =•  1 68. 570  (  pcf  ) 

I 

THE  ELASTIC  MODULUS  OF  SOLID  FOAM  =  £35.000  (  Lei  ) 

THE  YIELD  STRENGTH  OF  SOLID  FOAM  =  18.415  <  Lsi  ) 

THE  MASS  DENSITY  OF  SOLID  FOAM  =  75.000  C  pcf  ) 

THE  PROPORTIONAL  IT  Y  CONSTANT 

FDR  ELASTIC  MODULUS  OF  FOAM  CORE  =  1.130 

THE  POUER  CONSTAHl 

FOR  El  AST  1C  MODU!  US  OF  FOAM  CORE  =  1.710 

THE  PRUF'ORT  lONAL]  1  Y  CONSTANT 

FOR  SFIEAR  LIODUl  IJS  Dr  F(iAM  Cl.iRE  =  .400 

THE  POWER  CONSTANT 

FOR  SEAR  MLiDUl  US  OF  FOAM  CORE  =  £ ,  OO'j 

THE  RRDPDRTIDNALl'Ry  CONSTANT 

■^DR  SHEAR  STRENGTH  OF  FOAM  CORE  =  .310 

THE  POWER  CONSTANT 

FOR  SHEAR  STRENGTH  OF  FOAM  CORE  =  1.520 

SIMPLY-SUPPORTED  UNDER  CONCENTRATED  LUAr) 

THE  WIDTH  OF  THE  SANDWICH  BEAM  =  £.000  (  in.  ) 

THE  SPAN  OF  THE  SANDWICH  BEAM  ==  30.000  (  in.  ) 

THE  REOUIRED  DESIGN  STRENGTH 

PER  UNIT  WIDTH  A1 ’I  LENGTH  =  m.OOum  r;.Ej 

THE  DPT  I  HAL  DESIGN  v.'Al.U£S  ; 

THE  OPTIMAL  CORE  THICKNESS  = 

THE  optimal  r«CE  THICKNESS  ^ 

THE  OPTIMAL  NhSS  DENSITY  OF  FOAM  = 

THE  MINIMUM  WEIGHT  OF  SANDWICH  DEAN 


2.6S9E  in. 

.'0250  in. 
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.  .C"006  lopB 


FY/FW  FAILURE  DESIGN 


r'lINIMi.ll'l  WEIGHT  DESIGW  FOR  STIFFNESS  IH  SHiaDWlCH  BEAriS 
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MINinUM  UEIGHT  DEblGU  FOR  STIFFI-IESB  A  !f.  STREI'JGTH  IN  SANDWICH  EfEANS 
THE  ELASTIC  MODULUS  DF  FACE  MATERllAL  10150.000  <  Lei  ) 

THE  YIELD  STRENGTH  OF  FACE  MATERIAL  =  1E.A70  (  ksi  ) 

THE  MAES  DENSITY  OF  FACE  MATERIAL  =  16S.570  (  pcf  ) 

i 

THE  ElABTIC  MODULUS  OF  SOLID  FOAM  =  £32.000  (  Tsi  ) 

THE  YIELD  ST  RETTGiH  OF  SOLID  FC/Ari  =  1S.A15  <  Le:  ) 

THE  MASS  density  OF  SOL  1 0  FOAM  =  75.000  (  pcf  ) 

T  HE  RROF'UF:!  1 OMAL  ]  T  V  CONST  ANT 

FOR  ELASTIC  MDOUl  OS  OF  FOAM  CORE  =-  1.130 

tt-;e:  power  cdnstani 

FOR  elastic.  MGT'.Ui  I..L-  OF  FOfM  CORE  =-  1.710 

THE  PROPORT]  DNAI  j  1  y  CONST  FiNT 

FOR  SHEAR  MODULUS;  OF  FOAM  CORE  =  .AO'.i 

THE  POWER  CONE'.ThITI 

FOR  SEAR  MODULUS  OF  Fi  !(.•-,  n  Ci  EET  =  2.000 

THE  FRDRDRTIONAL] TV  CONS TAUT  '  ■ 

FOR  SHEAR  STRENGTH  DF  FOAM  CORE  =  .310 

THE  ROWER  CONSTAm 

FOR  SHEAR  STRENGTH  OF  FOAM  CORE  =  1.520 

SIMPLY- SUPPORT  ED  UNDER  CONCENTRATED  i  C  Ai' 


THE 

WIDTH  OF 

THE  SANDWICH  BEAM  = 

£ .  000 

( 

i  n 

THE 

SPAN  OF 

THE  SANDWICH  BEAM  = 

30 . 000 

< 

i  n 

THE 

REQUIRE)) 

DESIGN  LOAji  = 

.  E50T) 

1- 

■*-  P-' 

THE 

REQUIRED 

DESIGN  DEFL..ECTION  = 

.  50)00' 

3. 

n . 

THE  OPTIMAL  DESIGN  VALUE  ; 

THE  OPTIMAL  CORE  THICKNESS  = 

THE  OPTIMAL  FACE  THICKNESS  = 
THE--OPTIMAL  MASS  DENSITY  OF  FOAM  = 

THE  MINIMUM  WEIGHT  DF  SANDWICH  IEaVT  = 
E  i  IFFNESS.'F  >  -AILURE  DESIGN 


£.9965  in. 

. 0E5 1  in. 
6.5763  pcf 
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At4ALYSIS  OF  SAI-JDWICH  PLATES 

THE  ELASTIC  MODULUS  OF  FACE  MATERIAL  •  iO^Oo.OOO  (  i;5i  ) 

THE  YIELD  STRENGTH  OF  FACE  MATERIAL  =  IE. A 70  (  Lei  ) 

THE  MASS  DENSITY  OF  FACE  MATERIAL  =  '  16S.570  <  pcf  ) 

THE  ELASTIC  MODULUS  OF  SOLID  FOAM  =  E3E.000  (  Lei  ) 

THE  YIELD  STRENGTH  OF  SOLID  FOAM  =  18. A 15  (  Lei  ) 

THE  MASS  DENSITY  OF  SOLID  FOAM  =  75.000  (  pcf  ) 

THE  PROFORT]  ONAL  3  T  V  COOS':  A! 'T 

FDR  ELASTIC  Mui..U,„  ■  LL  OF  Fi  <Ai  1  COFC  ^  1.3  3':' 
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RECTANGULAR  PLAl E .  SIMPLY-SUPPORTED  UNDER  UNIFORM  LOAD 

THE  POISSON'S  RATIO  OF  THE  FACE  MATEF.IAL  .300 

THE  LENGTH  (  X  DI -:ECT  JON  >  OF 

THE  REC I  riNGMi.,  mK  F'LATE  =  Sl'0 . 0O'.>  (  iv-|.  ) 

THE  .WIDTH  (  V  DIRECTION  )  OF 

THE  RECTANGULAR  PLATE  =  100.000  (  in.  ) 

THE  MAGNITUDE  OF  THE  UNIFORM  LOAD  =  .00100  (  Lsi  ) 

THE  THICKNESS  OF  FACE  MATERIAL  =  .100  (  in.  ) 

THE  THICKNESS  OF  FOAM  CORE  =  4.000  (  in.  ) 

THE  MASS'^DENSTT’,  0-  FOAM  CORE  =  .  4.000  (  pcf  ) 

THE  MAXIMUM  STRESSES  IN  THE  FACE  MATERNAL  : 


THE 

MAXIMUM 

STRESS  BIGMAX 

( S I GMAR  > 

1 . 1300 

L,=  i 

THE 

MAX  I  liUM 

STRESS  SIBMAY 

(SIGMhT ) 

8. 4798 

t'Ei 

THE  MAXIMUM  STRESS  TAUXY  (TAURT) 


1  .  186^+ 
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THE  MAX  1  HUM  DEFI.ELTIDN  = 

THE  FAILLIF.'E:  MDOE  'S'  CORE  SHEA-Fi  YiELER'lE 
I  HE  SANDWICH  Pl.AiTE  IJNDEF  LOADINUi  IS  S-  FE 


K  '  ^  tZ' 

.  01  1  1 


1  .  1  t'OO 


i  E  1 
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MltJinilM  WEIGHT  rESIGH  T..  .  STIFFNESS  IN  SANDWICH  PLATES 
THE  ELASTIC  NDDULUS  L'F  FACE  IL,TE:F;IAL  lOlCO.OO'j  ’  i  zi  ) 

THE  YIELD  STRENGTH  OF  FACE  MATERIAL  =  1S.A70  L  =  i  ) 

THE  MASS  DENSITY  OF  FACE  MATERIAL  =  '  16&.L70  <  pcf  ) 

THE  ELASTIC  MODULUS  OF  SOi-ID  FOAM  =  S3R.000  l  ksi  ) 

THE  YIELD  STRENGTH  OF  SOLID  FOAM  =  IB. 415  (  ksi  ) 

THE  MASS  DENSITY  OF  SOLID  FOAM  =  75.000  (  pet  > 

THE  PROPORT IDNALITY  CONSTANT 

FOR  ELASTIC  MODULUS  OF  FOAI'l  CORE  =  1.130 

THE  POWER  CONSTf'TiT 

FOR'  EI-AST  IC  fT(IiL/ULU;6  Ur  FOAr'i  CORE  =  1 .71'.i 

THE  PROPORTIONALITY  CONST AMT 

FOR  SHEAR  MODULUS  OF  FOAM  CORE  =  .400 

THE  POITER  CONST  AM  i 

FOR  SEAR  MODULUS  OF  FC*AM  CORE  =  ’  E.OOO 

THE  FROPORT 1 ONAL. ;  T  V  CONSTANT 

FOF<  SHEAR  STRENGTH  OF  FOnM  CORE  =  .310 

THE  POWER  CQNSTAriT 

FOR  SHEAR  STRENGTH  OF  FOAM  CORE  =  ■  1.5S0 

RECTANGULAR  PLATE,  SIMPLY-SUPPORTED  UNDER  UNIFORM  LOAD 

THE  POISSOFL'S  RATIO  OF  THE  FACE  MATERIAL  .300 

THE  LENGTH  <  X  DIF'ECTION  )  UF 

THE  R.'E CT At'lGUi.  riF!  PL. ATE  =  200.000  '  i  ri .  ) 

THE  WIDTH  <  Y  DIRECTION  )  OF 

THE  RECTANGULAF;  PLATE  =  100.000  (  in.  ) 

THE  REC'UIRED  DESIGN  STIFFNESS  =  .00100  <  ksi/in  ) 

THE  OPTIMAL.  DESIGri  VALUES  ; 

THE"'oPTIMAI  THICl:NEBS  of  FACE  material  =  .0178  ill. 

THE  OPTIMAL  THICKNESS  OF  CORE  MATERIAL  =  5.53-'^'h  in. 

THE  OPTIMAL  MASS  DENSITY  OF  FOAM  =  4.5883  pcf 


THE  MINIMUM  WEIGHT  OF  SANDWICH  PLATE 


3464 


k:  1  p  E 


MIMIMUK  WEIGHT  DESIGT!  FOR  STRENGTH  IN  SArF-.HICH  PLATES 


THE  ELASTIC  MODULUS  C.R  FACE  MATERIAL  = 
THE  YIEl  D  STRENGTH  GF  FACE  MATERIAL.  = 
THE  MASS  DENSITY  OF  FACE  MATERIAL  = 

THE  ELASTIC  MODULUS  OF  SOLID  FOAM  = 

THE  YIELD  STRENGTH  GF  SOLID  FOAM  = 

THE  MASS  DENSITY  OF  SOLID  FOAM  = 


0150 

..  OC.i'O 

<  !:£:i 

) 

IH 

.  4  7  '0 

(  Lei 

) 

> 

168 

.570 

<  pcf 

) 

S3E 

.  000 

(  Lei 

) 

18 

.415 

(  L&i 

■) 

75 

.  000 

(  pcf 

) 

THE  PROPORT 1 DNAL I  TY  C GtTST  ALL! 


FDR  ELASTIC  MODULUS  C"  FOAM  CORE  = 

1  . 1  St* 

THE  POWER  CONSTANT 

FOR  El;  AS  TIC  MODULUS  0-'  FOAI';  CORE  = 

1  .710 

T!  IE  RROPQRT'  1  DNAL  J  T  V  C  CiNST  ANT 

FDR  SHEAR  Mi  iDi  )L LIB  OF  P  DAM  COT-E  = 

.  400 

THE  POWER  CONSTANT 

FOR  SEAR  MODULUS  0-  FuAM  CORE  = 

2 . 000 

THE  PROPORTIONALITY  CONSTANT 

FOR  SHEAR  STRENGTH  OF  FOAM  CORE  = 

.310 

THE  ROWER  CONSTANT 

FOR  SHEAR  STRENGTH  OF  FOAM  CORE  = 

1 . 520 

RECTANGULAR  PLATE,  S3 MPLY -SUPPORT ED  UNDEF 

UNIFORM  LOAD 

THE  POISSON'S  RATIO  OF  THE  FACE  MATERIAL 

.  300 

THE  LENGTH  (  X  DIRECTION  )  OF 

I  He  RE(JTANi->ULAr!  F'l,.ATE  = 

200.00CI  \  111 

.  ) 

THE  WIDTH  (  Y  DIRECTION  >  OF 

THE  RECTANGULAR  PLATE  = 

1  C>0 . 000  (  i  n 

.  > 

THE  REQUIRED  UNIFORM  LOADING  = 

1.0000  psi 

THE  OPTIMA!  DESIGN  VALUES  : 

THE  OPTIMAL. CORE  THICKNESS  = 

2 . 0729 

in. 

THE  OPTIMAL  FACE  THICKNESS  = 

.0179 

1  ri . 

THE  OPTIMAL  MASS  DENSITY  OF  FOAM  = 

2.^150 

pcf 

THE  MINIMUM  WEIGHT  OF  SANDWICH  PLATE 

=  .  1  :-)99 

1:  2  0  = 

FY/FW  FAILURE  DESIGN 


tIIMJriUM  UEIGHT  DEL  r .[  b'-l  FDn:  STIFFNESS  r-,  .  STRENGTH  IN  BANiLW  I CH  PLATES 
THE  ELASTIC  IIODULtS;  DF  FACE  riATERlAL  10150.000  <  Lei  ) 

THE  YIEl.D  STRENGTH'  D-  FACE  K,:,T ER I AL  1S.A70  (  l.si  > 

THE  NASS  UENSn  v'  tO’  FACE  MATERIAL  1  EE!.  570  (  pcf  ) 

THE  ELASIIC  NQDULl.'S  OF  EDLi)t  FOAM  =  S32.000  (  F:si  ) 

THE  YlEl.O  STRENGTH  OF  SOLID  FDAN  =  IB.-hIS  (  Lej  ) 

THE  MASS  DENSITY  OF  SOLID  FOAM  =  75.000  (  pcf  ) 

THE  FRDFORTinNAL.n  V  COriS  TArTI 

FOR  ELASTIC  MODULL'S  OF  FO-,M  CORE  =  1.130 

T  F^E  F'DWER  CCjI'IL' I  A! .!  i 


F  riR 

Ei  AS TIC;  MJOUL  L'.' 

OF  FL 

.'AM  L.Y.lr-;L  - 

1.71 

T  1  LS 

RRGPOR  1  1  UITAL  1  T 

CCNFMrAMT 

FOR 

SHEAR  MUDUi  US  (j-- 

F  UAI' 

i  COF-.'L  -- 

.  H  *  jO 

THE 

F'OR 

PUWLR  COIYSTAMT 
SEAR  MODULUS  OF 

FCiAr-: 

CORE  == 

E .  OCiiT 

T HE  ■ ROFORT I GNAL 1 7  V  CONSTANT 

FDR  SHEAR  SI RENGTH  OF  FDAN  CORE  =  .310 

THE  POWER  CONSTAfn 

FOR  SHEAR  STRENGTH  OF  FOAM  CORE  =  1 . 5E0 

RECTANGUl.AR  PLATE,  SIMPLY-SUPPORTED  UNIER  UNIFORM  LOAD 

THE  POISSON'S  RATKD  OF  THE  FACE  MATERIA!  .3  00 

THE  LENGTH  (  X  DIRECTION  )  OF 

THE  RECTANGULAR  PLATE  =  EOO.OOO  (  ■ 

THE  WIDTH  (  Y  DIRC'  'HON  ’  OF 

T  F'3  Rl-.CT  IL '  ■  F'L.  i  £  =  I  '.lu .  (  i  ri .  ) 

THE  REQUIRED  UNIFORM  LOADING  =  .0010  Lei 

THE  DESIGN  DEFLECTION  =  1.0000  in. 

THE  OPTIMAL  DESIGN  VALUE  : 

THE  OPTIMAL  CORE  THICKNESS  = 

* 

THE  OPTIMAL  FACE  THICKNESS  = 

THE  OPTIMAL  MA3 S  DENSITY  OF  FOAM  = 

THE  MINIMUM  WiE  I GHT  OF  SANDUIOH  PLAT  i 
STIFFNESS/PV  F^■,:LL. IJRE  DESIGN 


1 C) .  1 E07  1  n . 

. 0037  i n . 
3.o316  pc.f 
.^397  i:-jpE 
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THICK^ESS  OF  FOAM  CORE 


DENS  I T 


